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ABSTRACT 


The  results  of  an  expanded  study  and  investigation  of  the 
Pulse  Pair  technique  for  estimating  the  first  and  second 
moments  (mean  and  variance)  of  doppler  spectra  for  radar 
backscatter  from  atmospheric  phenomena  are  presented. 
The  theory  is  extended  tc  include  the  effects  of  non-ideal 
conditions,  such  as  noise,  and  experimentally  verified  by 
extensive  performance  tests  using  simulated  weather 
signals  with  controllable  parameters.  A  proposed  experi¬ 
mental  model  of  a  real-time  digital  pulse  pair  processor 
is  defined  and  compared  with  alternate  processing  tech¬ 
niques.  Based  on  the  encouraging  results  of  the  study, 
recommendations  are  made  to  carry  the  theory  into 
practice;  these  include  the  construction  of  a  real-time 
digital  pulse  pair  processor  with  flexible  characteristics 
to  gather  and  reduce  data  for  evaluation  while  operating 
with  real  radars,  and  the  development  of  additional 
related  theory  needed  to  guide  the  experimental  effort. 
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1.0  INTRODUCTION 


1.  I  Scope 

This  la  the  Final  Report  on  Contract  F  19628-71-C-0126,  which  is 
devoted  to  the  study  and  investigation  of  the  pulse  pair  measurement  tech- 
niqu  for  determination  of  the  first  and  second  moments  of  Doppler  spectra 
of  radar  backscatler  from  atmospheric  phenomena.  This  technique  pro¬ 
mises  to  greatly  simplify  and  thereby  reduce  the  costs  of  deriving  quantita¬ 
tive  data  from  a  pulse  Doppler  weather  radar.  Initial  efforts  in  the  applica¬ 
tion  of  the  pulse  pair  technique  to  weather  radar  signals  were  conducted 
under  a  portion  of  a  prior  contract  with  the  Weather  Radar  Branch  of  the 
AFCRL  Meteorology  Laboratory  (F 19628-68-C-0345).  The  emphasis  on  the 
present  contract  has  been  to  extend  the  theory  and  to  experimentally  verify 
the  performance  of  the  pulse  pair  estimating  technique  to  include  the  effects 
of  non-ideal  conditions,  such  as  noise,  that  are  encountered  in  an  operational 
environment  with  real  equipment;  and  to  prepare  a  design  specification  for  a 
real  time  pulse  pair  signal  processor  that  could  be  utilized  in  weather  radar 
research. 

In  the  interests  of  overall  continuity  and  perspective,  and  to  facilitate 
understanding,  this  report  provides  a  unified,  self-sufficient  presentation  of 
the  theory  and  experimental  work  to  date,  including  some  of  the  early  results 
from  the  previous  contract. 

1 . 2  Summary 

Major  steps  forward  *hat  were  taken  during  the  present  contract  include 
(1)  establishing  a  theoretical  foundation  for  quantitatively  predicting  the  per¬ 
formance  capabilities  of  the  pulBe  pair  estimators  under  various  combinations 
of  conditions,  and  optimizing  the  design  parameters  of  the  radar-processor 
system;  and  (2)  experimentally  verifying  the  theory  by  hundreds  of  quantitatively 
controlled  performance  tests  conducted  with  the  aid  of  the  simulation  system 
developed  for  this  purpose. 

These  studies  have  demonstrated  a nalytically,  the  feasibility  of  im- 
ple  menting  a  real  time  signal  processor,  using  the  pulse  pair  computing 
algorithm  as  its  central  core,  to  obtain  quantitative  estimates  of  mean  velocity 
and  velocity  spread  under  a  wide  variety  of  conditions  and  for  a  wide  variety 
of  radars.  This  report  describes  the  initial  concepts  and  features  of  a  processor 
enabling  experimental  confirmation  of  its  predicted  performance,  in  a  form 
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suitable  for  performing  meteorological  research  with  Doppler  radars  as  well. 
With  the  instrument  proposed  here  it  will  be  possible  to  collect  and  reduce 
data  under  a  wide  variety  of  conditions  economically.  This  will  make  it 
possible  to  establish  a  broad  data  base  from  which  the  operational  utility  of 
a  Doppler  weather  radar  may  be  inferred.  The  proposed  signal  processor 
will  then  provide  a  benchmark  for  sizing  and  designing  a  signal  processor 
for  operational  use. 

The  key  features  available  in  a  pulse  pair  signal  processor  are 

1.  quantitative  (digital)  estimates  of  spectral  mean,  spectral  width 
and  target  reflectivity  in  real  time 

2.  large  numbers  of  contiguous  range  gates  may  be  processed 
simultaneously  (up  to  1024  cells  in  the  proposed  implementation) 

3.  selectable  dwell  time  may  be  used  to  adjust  the  block  length  of  the 
signal  processor  to  the  radar  scanning  requirement  and  to  the 
meteorological  conditions  (integration  of  up  to  1024  samples  per 
range  cell  to  be  provided  in  proposed  signal  processor) 

4.  range  resolution  may  be  obtained  (resolution  equivalent  to 
ft  (.  5 ijis )  in  the  proposed  implementation) 

5.  predictable  performance  based  on  self  contained  data  (confidence 
bounds  on  all  estimates  are  computable  based  on  signal  processor 
derived  S/N  data) 

6.  sliding  window  operation  is  feasible  (estimates  based  on  the 
most  recent  N  seconds  of  data) 

7.  the  storage  required  consists  essentially  of  three  complex  words 
per  range  cell  (last  sample,  summary  crosscorrelation  data, 
summary  power  data)  while  the  computing  power  required  is  to  be 
able  to  calculate  2  complex  multiplies  and  adds  at  the  input  data 
rate 

8.  additional  waveform  flexibility  (non  uniform  p.  r.  f.  and  pulse  -'idths 
limited  frequency  agility 

9.  sensitivity  comparable  to  spectrum  analyzers 

10.  independent  AGC  for  each  range  bin  to  improve  overall  dynamic 
range 
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as  well  as  the  features  common  to  all  digital  signal  processors,  namely 

11.  non  critical  adjustments 

12.  long  term,drift  free, stable  performance  without  calibration 

13.  flexibility  in  interfacing  with  a  variety  of  radars. 

The  digital  outputs  of  the  processor  can  be  displayed,  recorded,  or 
remoted  for  further  display  or  processing. 

To  fully  utilize  the  processor  capabilities  it  must  operate  in  conjunction 
with  a  coherent  radar  receiver  that  has  a  wide  dynamic  range  and  is  able  to 
accept  AGC  control  voltages  from  the  processor  at  a  range  cell-to-range 
cell  rate.  To  insure  compatibility,  it  is  proposed  to  incorporate  in  the 
processor  a  self-contained  receiver  with  the  required  characteristics  that 
will  accept  1.  F.  and  reference  signals  from  the  radar,  gain  control  voltages 
from  the  processor,  and  deliver  the  required  in-phase  and  quadrature  complex 
video  signals  to  the  processor.  Since  the  processor  measures  signal  amplitudes 
to  form  the  AGC  control  voltage  for  each  range  cell,  this  information  is  avail¬ 
able  as  a  digital  output.  Together  with  the  large  useful  receiver  dynamic 
range,  as  an  added  benefit  the  processor  amplitude  output  in  conjunction  with 
range  information  can  also  be  utilized  for  other  purposes,  such  as  the  cal¬ 
culation  of  target  reflectivity. 

It  should  be  noted  that  although  there  are  other  methods  of  estimating 
spectral  mean,  the  pulse  pair  technique  provides  quantitative  estimates  of 
this  parameter  in  digital  form,  and  is  the  only  known  method  for  estimating 
spectral  width  in  real  time. 

The  balance  of  this  report  is  devoted  primarily  to  the  technical  aspects 
of  the  pulse  pair  program.  The  technical  background  is  presented  in  the 
succeeding  portion  of  this  Introduction;  the  theory  of  Pulse  Pair  Estimation, 
including  original  analyses  performed  on  the  present  contract,  appears  in 
Section  2;  Section  3  deals  with  the  methods  and  results  of  Experimental  Eval¬ 
uations  of  the  pulse  pair  technique,  including  that  performed  with  real  radar 
data  on  the  prior  contract  and  that  utilizing  simulated  radar  return  on  the 
present  contract;  Section  4  describes  the  proposed  Hardware  Realization 
of  a  Pulse  Pair  Processor  capable  of  operating  on  real  radar  data  in  real 
time;  and  Section  £  presents  the  conclusions  that  have  been  drawn  and  the 
Recommendations  that  have  oeen  formulated  as  a  result  of  the  Study  and 
Investigation. 
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Detailed  Analyses,  Computer  Programs,  and  Spectral  Plots  are  in¬ 
cluded  in  the  Appendices  in  support  of  the  main  text  of  this  report. 

1.  3  Technical  Background 

A  commonly  occurring  problem  in  radar  measurements  is  to  estimate 
the  spectral  mean  and  variance  of  a  unimodal  Gaussian  signal  imbedded  in 
white  Gaussian  noice.  The  classical  technique  used  is  to  form  spectral 
estimates  from  the  time  sequence,  and  then  form  estimates  of  the  spectral 
mean  and  variance  by  standard  formulae.  It  is  not  difficult  to  show  that  these 
techniques  are  not  optimum  and  result  in  signal  processor  dependent  bias 
errors  which  affect  the  system  sensitivity. 

When  both  mean  and  variance  are  unknown,  the  optimum  spectral  para¬ 
meter  estimates  can  only  be  obtained  implicitly  from  the  solution  of  a  trans¬ 
cendental  equation  in  which  the  form  of  the  underlying  spectrum  is  known. 

To  overcome  these  defects,  various  investigators  have  resorted  to  threshold 
tests  to  reject  the  noisy  frequencies,  in  order  to  apply  classical  mean  and 
variance  estimators  to  the  residue  after  thresholding. 

This  report  describes  a  new  technique,  known  as  pulse  pair  analysis, 
which  bypasses  the  stage  of  spectrum  analysis,  and  directly  estimates  spectral 
mean  and  variance.  Not  only  docs  this  technique  result  in  a  less  complex  set 
of  calculations  to  be  made  by  the  processor,  but  it  has  a  form  which  is  based  on 
an  optimum  estimation  theory  and  as  such  is  relatively  free  of  arbitrary  choices. 
Furthermore,  the  technique  is  robust  with  respect  to  underlying  spectral 
shape,  and  permits  more  flexibility  in  the  choice  of  radar  waveform. 

Its  major  disadvantage  with  respect  to  classical  spectral  parameter  esti¬ 
mation  is  the  fact  that  the  technique  does  not  indicate  the  presence  of  multi- 
peaked  spectra.  When  the  full  spectral  analysis  is  available  clutter  regions 
may  be  identified  and  possibly  rejected,  and  separated  non-standard  spectral 
shapes  may  be  detected.  Pulse  pair  analysis,  being  a  spectral  shape  inde¬ 
pendent  estimation  technique,  requires  a  clutter  rejection  device  preceding  it 
(where  necessary),  and  yields  only  indirect  evidence  of  peculiar  spectrum 
conditions  (abnormally  broad  widths)  for  indicating  the  presence  of  multiple 
lobed  spectra. 

As  for  sensitivity,  simulation  evidence,  to  be  presented  in  a  later  sec¬ 
tion,  indicates  that  the  sensitivity  of  the  pulse  pair  estimation  is  comparable 
to  that  of  the  shape  independent  spectrum  processing  techniques. 
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2,  0  THEORY  OF  PULSE  PAIR  ESTIMATION 


2.  1  The  Pulse  Pair  Estimators 

It  may  be  shown  ^  that  the  optimum  estimate  of  the  spectral  mean  and 
variance  cf  a  sequence  of  pairs  of  measurements  which  are  uncorrelated  is 
given  by 
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rik  is  the  observed  signal,  S(f)  is  the  underlying  spectral  density 
function,  and  is  the  mean  equivalent  noise  cross-section  (see  Figure  1). 

It  is  possible  to  rewrite  these  equations  in  a  more  illuminating  form  by 
defining 


We  get 


R(t)  =  y*  S(f)  e 


J2nft 
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Thus  the  latter  equation  alone  determines  w,  and  furthermore,  if  R(T)  is  real, 

A 

the  first  equation  is  independent  of  w. 

Essentially  this  result  says  that  if  the  underlying  spectral  density  function 
of  the  signal  is  characterized  by  a  spectrum  whose  shape  is  invariant  to  frequency 
shift  and  broadening  in  the  specified  manner,  then  these  two  parameters  may  be 
estimated  from  the  pair  wise  measurement  directly  by  estimating  the  complex 
correlation  function  of  the  data. 


It  is  important  to  note  that  for  small  T 

R(wT)  ?  R<0)-1/2wV'(0)  (2irT)2 


(2-8) 


(No  first  order  term  is  present  because  S(f)  by  definition  has  mean  zero. )  Also  by  defi¬ 
nition  R(0)  and  R"  (0)  are  both  unity  (the  first  to  have  unit  energy,  the  second 
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to  define  the  parameter  w,  i.  e. ,  unit  variance  for  the  underlying  signal  shape). 

As  a  result,  in  this  case,  the  estimator  of  spectral  width  parameter  is  given 

by  _ 
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and  the  estimator  of  the  spectral  mean  is  given  by 


*  2rTT  are  x‘ 


(2-10) 


These  are  the  pulse  pair  spectral  parameter  estimators  suggested  by 
Rummler.  (2)(3)(4) 


The  approximations  made  in  arriving  at  these  results  indicate  that  no 
difficulty  in  utilizing  these  estimators  should  arise  beyond  that  expected  from 
sampling  considerations.  Thus  frequencies  in  excess  of  will  be  ambiguously 
folded  by  the  mean  estimator,  and  the  variance  estimator  will  be  degraded  when 
the  underlying  spectral  width  begins  to  approach  the  sampling  frequency. 
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It  is  also  important  to  note  that  the  explicit  estimators  do  not  depend 
directly  on  any  particular  spectral  shape.  This  is  not  true  of  estimators  of 
these  same  parameters  based  on  spectral  measurements. 


2.  2  Estimation  Accuracy 

Since  (2-9)  and  (2-10)  are  explicit  functions  of  the  data,  it  is  straightforward 
to  analyze  the  expected  performance  of  the  estimators.  The  analysis  is  com¬ 
plicated  by  the  nonlinearity  of  the  estimator  functions. 


Appendix  A  gives  an  approximate  analysis  of  the  performance  of  both  the 
mean  and  the  width  estimators.  It  is  shown  there  that  the  estimator  is  unbiased 
and  that  when  the  per  puloe  signal  to  noise  ratio  is  useful,  the  variance  of  these 
estimators  is  given  by 


2M  Var  (2irfQT)  =  e(2lTwT) 
and 

2M  (2ttwT)2  Var  (2ttwT)  = 


+  2irwT 
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where 

M  =  number  of  pairs  in  the  estimate  («  number  of  data  points  in  the 
sample) 

N/S  =  per  pulse  noise  to  signal  ratio 


These  formulae  were  derived  for  Gaussian  random  variables  with  Gaussian 
shaped  spectra.  More  general  formulae  are  available  in  Appendix  A. 


The  general  formulae  were  derived  for  a  uniform  train  of  samples  and 
account  for  all  correlations  in  the  data. 

(4) 

Figures  la  and  2  show  plots  of  these  functions.  Comparable  figures 
derived  for  the  case  of  independent  pulse  pair  measurements,  (i.e.,  from 
a  nonuniform  pulse  train,  with  successive  pairs  spaced  so  far  apart  that  the 
echoes  are  completely  decorrelated)  are  shown  in  Figures  3  and  4.  Comparison 
of  these  figures,  which  show  negligible  differences  at  useful  signal  to  noise 
ratios,  are  indication  again  of  the  robustness  of  the  technique. 
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Figure  2.  Standard  Deviation  of  Pulse  Pairs 
Width  Estimate  (wT)  va  True  Spectral  Width 
(wT)  with  Per  Pulse  S/N  a  Parameter,  for 
Uniform  Pulse  Spacing 


Spectral  Width  Parameter,  w 


re  4.  Fractional  Accuracy  of  Pulse  Pairs  Estimate  of  Spectral  Width 
True  Spectral  Width  for  Independent  Pairs  (from  Rummler,  Ref.  (4)) 


The  pulse  pair  technique  calculation  In  essence  requires  2  complex  pair 
multiplies  and  adds  for  each  new  data  point  in  order  to  form  X  and  Y  as  indicated 
in  (2-3)  and  (2-4).  The  conversion  of  these  to  mean  and  variance  estimates  is 
an  operation  which.need  not  be  carried  out  at  the  input  data  rate.  To  bring  this 
into  sharp  focus,  a  digital  spectrum  analyzer  capable  of  processing  up  to  1024 
complex  samples  at  a  single  range,  would  require  5  times  the  number  of  com¬ 
parable  arithmetic  operations  and  up  to  500  times  the  data  storage  of  a  com¬ 
parable  pulse  pair  processor. 


This  efficiency  in  signal  processing  makes  it  possible  to  consider  im¬ 
plementing  a  real  time  signal  processor  which  can  obtain  spectral  mean  and 
variance  estimates  in  real  time  with  selectable  dwell  time  per  gate  of  up  to  a 
1000  complex  samples  per  dwell,  with  as  many  as  1000  range  cells  processed 
simultaneously.  Available  technology  currently  permits  these  calculations 
to  be  made  at  a  2  MHz  data  rate,  thereby  permitting  range  resolution  0.  5^s 
(i.  e.  ,  250  ft.  )  and  even  higher  thruput  rate  processors  are  being  developed 
now. 

In  a  later  section  we  describe  a  candidate  signal  processor  implementing 
this  technique  to  the  above  performance  specifications.  That  section  details 
the  flexibility  and  performance  which  may  be  expected  of  a  relatively  modest 
signal  processor. 


2.4  Waveform  Flexibility 


The  restriction  of  the  processor  to  consecutive  pairs  makes  possible 
a  great  deal  of  waveform  flexibility  which  can  be  used  to  advantage  to 
overcome  certain  inherent  radar  limitations.  Prime  among  these  are  the 
range-Doppler  coupling  which  makes  the  unambiguous  Doppler  interval  of  the 
radar  vary  inversely  with  its  unambiguous  range  interval.  For  weather 


radars,  this  amounts  to  a  reduction  of  the  Doppler  analysis  band  of  the  radar 
as  its  range  is  extended. 


The  pulse  pair  technique  makes  it  possible  to  use  the  waveform  shown  in 
Figure  5a  to  change  the  range-Doppler  coupling  relationship  to 


(2-13) 
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a)  Simple  Pair  Spaced  Waveform 
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Figure  5*  Pulse  Pa.ir  Waveforms 


This  makes  it  possible  to  increase  the  Doppler  coverage  band  on  relatively 
small  colls  located  at  greater  distances  from  the  radar.  Selection  of 
and  permits  the  radar  to  operate  with  an  unambiguous  Doppler  spread 
reciprocally  related  to  the  dimension  of  an  isolated  cell.  In  weather  radars 
such  a  technique  might  permit  the  unambiguous  Doppler  measurement 
of  such  conditions  as  severe  convective  cells  and  possibly  local  CAT  con¬ 
ditions  at  low  angles  and  long  ranges. 

The  second  waveform  possibility,  shown  in  Figure  5b,  is  available  with 
both  equally  and  unequally  spaced  pairs.  The  accuracy  equations,  (2-11)  and 
(2-12),  reveal  that  a  key  performance  parameter  is  the  per  pulse  signal  to 
noise  ratio.  The  equations  indicate  that  there  is  an  optimum  S/N  per  pulse 
which  minimizes  the  measurement  errors  for  a  fixed  radar  energy.  That  is, 
if  S/N  and  M  are  related  by  the  equation 


£  =  M  (S/N)  (2-14) 

(which  indicates  that  a  fixed  energy  is  available  to  the  measurement  that 
may  be  distributed  in  any  way  between  the  number  of  pairs  processed  and  the 
S/N  used  per  pulse)  then  either  equation  (2-  1 1)  or  (2- 12)  has  a  minimum  at  a 
definite  value  of  M. 


In  particular,  it  may  be  shown  that  for  the  measurement  of  mean 
velocity,  the  optimum  per  pulse  S/N  is  given  by 

(5/Nlop,  =  12ir  wT  y[^)'1/2 

and  at  that  point  the  measurement  variance  is  given  by 
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Figures  6  and  7  indicate  the  sensitivity  of  this  condition  to  nonoptimum 
conditions  for  mean  Doppler  and  width  measurements  respectively.  It  is 
important  to  note  that  minimum  error  for  both  measurements  does  not  occur 
at  the  same  value  of  signal  to  noise  ratio.  Indeed,  the  curves  show  that  sub¬ 
stantially  more  per  pulse  S/N  is  required  to  minimize  width  errors. 


2-10 


l^mw-rr}  Va 


2-11 


idlduiuilJUiiKL  mlL>.hj.>.iK  ..uT,  •  .mii..  '  MallllilujJlii jatfaft  j*  liiwlLillilll 


The  tradeoff  characteristics  may  be  used  to  establish  the  benefits  to 
be  obtained  by  the  frequency  modulation  waveform  of  Figure  5b  as  follows. 

Suppose  that  the  radar  and  meteorological  parameters  are  such  that  only  normalized 
widths  of  0.05  or  greater  are  of  interest.  Then  Figure  7  indicates  that  the  data 
should  be  collected  at  the  signal  to  noise  ratio  of  7.  5  dB  per  pulse.  Suppose 
then  that  it  were  required  to  make  that  measurement  to  an  accuracy  of  10%. 

Then 
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For  wT  =  0.  05,  Figure  7  indicates  that 
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so  that 
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Thus, at  least  M  =  (36-7.5)dB 
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accuracy. 


=  930  pulses  must  be  processed  to  achieve  the  requisite 


If  the  peak  power  in  the  pulse  were  such  that  13.  5  dB  per  pulse  were 
available  for  the  measurement,  this  theory  suggests  that  transmission  of 
four  frequency  coded  subpulses  each  with  /.  5  dB  per  pulse  signal  to  noise 
ratio  for  a  total  transmission  of  250  groups  of  four  might  result  in  bette 
overall  accuracy  than  would  the  same  length  of  uncoded  transmission.  The 
improvement  in  performance  would  be  negligible  at  wT  =  .05  for  the  con¬ 
ditions  given,  but  at  wT  =  0.  1  the  frequency  modulation  would  improve  the 
accuracy  by  40%  ts  (■yj  2-1). 
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There  is  an  important  caveat  which  must  be  considered  in  this  regard; 
namely,  that  additional  receivers  may  be  required  to  implement  the  frequency 
modulation  system  for  extended  targets. 

2.  5  Summary 

The  pulse  pair  technology  offers  new  potential  for  Doppler  weather  radar 
systems.  It  permits  the  economic  implementation  of  a  real  time  signal  processor 
which  can  produce  prodigious  data  reduction,  reducing  the  enormous  amount 
of  data  produced  by  such  a  radar  to  more  manageable  dimensions.  It  does  so 
without  a  sacrifice  in  sensitivity  in  cases  of  interest. 

It  makes  it  possible  to  consider  new  radar  modulation  techniques  per¬ 
mitting  the  radar  to  collect  the  data  more  efficiently.  Thus  a  frequency  and  pulse 
pair  spacing  agile  radar, which  has  merit  in  both  increasing  the  data  rate 
(number  of  pulse  pairs  available)  and  suppressing  some  of  the  range-Doppler 
ambiguities, can  be  processed  without  penalty  in  a  pulse  pair  analyzer. 

The  pulse  pair  technique  produces  estimates  which  have  predictable 
accuracy  and  are  free  of  arbitrary  parameter  selection  (such  as  the  threshold 
level  in  the  spectrum  analyzer  technique).  In  a  digital  processor,  the  block 
length  of  the  analysis  interval  may  be  varied  on  command  so  that  the  processor 
can  easily  be  properly  matched  to  the  radar  operating  conditions  (p.  r.  f.  and 
scan  rate). 

If  data  rates  were  not  a  consideration,  it  is  clear  from  the  structure  of  the 
processor,  that  it  could  be  implemented  as  a  sliding  window  processor,  so  that 
running  averages  of  the  most  recent  N1  pulse  pair  data  could  be  used  to  obtain 
continuously  updated  mean  and  variance  estimates. 

In  short,  it  appears  that  the  pulse  pair  processing  technique  offers  an 
efficient  means  of  automating  pulse  Doppler  weather  radar  processing  with  a 
technique  which  is  accurate,  sensitive  and  flexible.  The  equipment  used  in 
conjunction  with  clutter  cancelling  techniques  should  permit  operation  in  a 
significant  background  of  ground  clutter. 

Its  ability  to  signal  the  presence  of  unusual  spectral  shape  is  perhaps 
its  ultimate  limitation. 
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3.0  EXPERIMENTAL  EVALUATION 


This  section  describes  several  efforts  to  date  to  experimentally  verify 
the  theoretically  predicted  accuracy  capabilities  of  the  pulse  pair  estimators 
of  spectral  mean  and  spread  and  to  compare  their  performance  with  other 
methods  of  estimation  using  identical  input  data.  The  initial  experiment,  which 
was  conducted  on  the  previous  contract,  compared  the  results  of  pulse  pair  and 
conventional  spectral  analysis  processing  of  real  Porcupine  Radar  weather 
return^'  The  limitations  encountered  with  this  approach  to  evaluation  moti¬ 
vated  the  more  extensive  and  precise  experimental  evaluations  accomplished  on 
the  present  contract  using  simulated  radar  return.  Although  all  of  the  experi¬ 
ments  to  date  have  involved  non-real  time  pulse  pair  processing  of  the  tape 
recorded  radar  output  data,  real  time  processing  is  both  feasible  and  advan¬ 
tageous;  a  proposed  real  time  hardware  processor  is  described  in  Section  4. 

3.  1  Porcupine  Radar  Weather  Return 

A  sequence  of  10240  consecutive  raw  complex  video  samples  (both  in-phase 
and  quadrature)  were  tape  recorded  from  the  output  of  the  Porcupine  weather 
radar  during  a  light  rain.  The  PRF  was  3300  pps,  the  antenna  elevation  angle  was 
40“  ,  and  the  data  was  obtained  from  the  range  cell  with  the  maximum  S/N. 

The  analog  data  was  then  quantized  for  subsequent  processing  on  a  digital  com¬ 
puter  by  the  "pulse  pair"  method,  with  conventional  spectral  analysis  as  a 
control.  The  latter  was  required  since  the  true  spectral  parameters  of  the 
radar  doppler  return  from  the  rain  were  unknown. 

3.  1,  1  Estimates  of  Mean  and  Spread  by  Conventional  Spectral  Analysis 

The  recorded  data  was  processed  by  conventional  spectral  analysis  in 
three  different  block  lengths:  one  block  of  8192  complex  samples;  two  blocks 
of  4096  samples;  and  ten  blocks  of  1024  samples.  A  typical  plot  of  the  amplitude 
of  the  spectral  components  of  one  of  the  ten  blocks  of  1024  samples  (representing 
approximately  1  second  of  weather  data)  as  analyzed  by  a  digital  computer  pro¬ 
grammed  to  execute  a  Fast  Fourier  Transform  is  shown  in  Figure  8,  "Typical 
Dat-  Spectrum  for  a  Sequence  of  1024  Samples.  "  The  rag,  ~d  appearance  of  the 
spectrum  is  due  partly  to  noise  and  partly  to  the  distributee^  nature  of  the  weather 
target. 
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The  computer  was  programmed  as  follows  to  estimate  the  first  moment 
(fQ)  of  the  doppler  spectrum  and  the  second  moment  (w^)  about  £  . 

tH 

Letting  z(n)  =  x(n)  f  jy(n)  denote  the  n  sample  in  a  block  of  length  N, 
the  spectral  amplitudes 
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are  computed,  where  the  a(n)  are  60  dB  Dolph  Chebyshev  weight  coefficients 
for  a  block  of  length  N. 

The  unambiguous  frequency  interval  is  then  shifted  from  0<  k  <  N-l 
N  N 

to  -  —  +  1  <  m  <  j  -  l  by  the  transformation 


Z(m)  = 


P(m-l), 


1  <  m  <  N  /l 
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as  shown  in  Figure  (4). 

The  spectral  mean  and  width  couldbe  calculated  in  accordance  with 


,  _  PRF  - 

in  ~"N~  m  »  wbere  m  = 


-  _  T,  mZ(m) 
Z(m) 


a  _  PRF  |  ^(m  -  m)*"  Z  (m)  ) 
"  N  V  SZ(m)  ] 


(3-3) 


(3-4) 


However,  these  classical  algorithms  for  the  mean  and  width  of  the  distribution 
make  no  provision  for  the  effects  of  noise,  although  it  is  ev'dent  from  Figure  8 

that  portions  of  the  spectral  plots  are  predominantly  i _ se,  which  would  tend  to 

obscure  and  degrade  the  estimates  of  the  spectral  parameters.  Neither  is  there 
a  developed  theory  for  the  effects  of  noise  on  their  accuracy  for  spectra  derived 
data.  Nevertheless,  despite  the  lack  of  guidance  from  theory,  it  is  necessary 
to  employ  some  additional  (though  not  necessarily  optimum)  technique  together 
with  the  spectral  analysis  estimators  to  improve  their  accuracy  so  that  they 
may  serve  as  standards  of  comparison  for  the  pulse  pair  estimators.  A  common 
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method  is  to  select  *  threshold  level  referenced  to  the  noise  level  or  to  the 
peak  of  the  spectrum,  and  to  use  only  those  spectral  components  whose  ampli¬ 
tudes  exceed  the  threshold  in  the  computation  of  spectral  mean  and  width  as  a 
special  set|  which  falls  within  the  general  class  of  spectral  analysis  estimators. 
Such  a  preliminary  computational  step,  prior  to  the  above  calculation  of  the 
mean  and  width,  has  been  introduced  into  the  spectral  analysis  estimator  com¬ 
puter  program. 

Rather  than  restrict  the  results  of  the  spectral  analysis  estimates  to  a 

A  A 

single,  possibly  poor  choice  of  threshold,  the  estimates  for  f  and  w  were 
computed  for  each  block  of  data  over  a  range  of  assumed  noise  thresholds  from 
20  dB  to  30  dB  below  the  spectral  peak,  and  are  shown  in  Tables  1  and  2  respec- 

A 

tively.  It  will  be  noted  that  the  estimates  of  f  vary  with  threshold  level  within 
each  block  over  a  range  of  from  5.  5%  to  12.  9%  of  the  nominal  value,  depending 
upon  the  specific  block. 

The  spectral  analysis  estimates  of  w  are  even  more  sensitive  to  the 
value  of  threshold.  They  vary  with  threshold  level  within  each  block  over  a 
much  larger  range  of  from  51%  to  68%  depending  upon  the  specific  block.  The 
potential  problem  of  using  such  measurements  as  a  standard  of  comparison  for 
evaluating  the  performance  of  the  pulse  pair  estimator  is  noted  here  in  passing, 
and  will  be  discussed  further  later  on. 

3.  1.  2  Estimates  of  Mean  and  Spread  by  Pulse  Pair  Processing 

The  identical  Porcupine  radar  data  processed  by  spectral  analysis  were 
also  processed  by  the  pulse  pair  technique.  The  ten  blocks  of  digitized  se¬ 
quences  of  1024  samples  were  inputted  to  a  computer  programmed  to  imple¬ 
ment  the  pulse  pair  estimators  of  equations  (2.  9)  and  (2.  10).  These  results 
are  also  shown  in  Table  1  and  2  respectively.  The  estimate  of  w  includes  a 
correction  for  noise,  the  level  of  which  was  deduced  from  the  experimental 
data  by  assuming  that  the  peaked  portion  of  the  spectrum  was  composed  of 
signal  plus  noise,  while  the  flat  skirts  were  noise  alone. 

3.  1.  3  Comparison  of  the  Estimators 

Numerical  comparisons  of  the  pulse  pair  and  spectral  analysis  estimates 

A  . 

of  f  and  w  obtained  from  the  Porcupine  weather  return  data  may  be  made  from 
the  results  in  Tables  1  and  2. 
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These  tables  show  the  values  of  f  and  $  as  estimated  block  by  block 

o  o 

by  pulse  pairs  and  for  each  threshold  level  for  the  spectral  analysis  method. 
The  blocks  of  1024  samples  have  been  numbered  sequentially  to  identify  the 
data  from  which  the  estimates  were  derived.  The  different  estimates  for  a 
given  block  of  data  appear  on  a  horizontal  line  alongside  that  block. 

As  was  noted  above,  the  spectral  analysis  estimate  of  f  varies  with  the 
threshold,  T.  It  agrees  best  with  the  pulse  pair  estimate  of  i  when  T  =  20  dB, 
for  which  case  the  agreement  is  within  better  than  1%  when  the  block  length, 

N,  is  4096  or  greater.  Since  the  spectral  analysis  estimates  of  fQ,  which  were 
used  as  the  control,  do  not  vary  too  widely,  it  was  generally  concluded  that 
the  pulse  pair  estimate  of  f  was  definitely  satisfactory  under  the  conditions 
for  which  the  data  was  taken. 

On  the  other  hand,  the  spectral  analysis  estimates  of  doppler  spread, 
w,  are  quite  sensitive  to  the  value  of  T.  Varying  as  they  do  over  a  range  of 
about  60%,  they  do  not  provide  a  satisfactory  standard  of  comparison  for 
evaluating  the  performance  of  the  pulse  pair  estimates  of  w.  Without  a  know¬ 
ledge  of  the  true  doppler  spread  of  the  weather  radar  return,  it  was  impossible 
to  draw  any  conclusions  as  to  the  effectiveness  of  the  pulse  pair  method  for 
estimating  with  this  particular  data. 

It  appeared  clear  from  the  above  results,  that  real  radar  weather 
return  of  opportunity  did  not  provide  a  promising  source  of  data  from  which 
the  behavior  of  the  pulse  pair  estimators  could  be  quantitatively  explored. 

Not  only  is  there  no  assurance  of  obtaining  data  with  the  range  of  spectral 
and  S/N  parameters  required  for  a  valid  exploration,  but  the  accuracy  with 
which  the  true  values  of  the  parameters  embedded  in  the  data  (which  are 
needed  as  a  basis  for  comparison)  can  be  established  is  in  itself  a  function 
of  the  parameters.  Since  we  cannot  reasonably  evaluate  one  unknown  against 
another  unknown,  the  alternate  approach  of  using  simulated  radar  return  with 
prescribed  parameters  was  indicated. 
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3.  2  Simulated  Radar  Return 

3.2.1  Purpose,  Capabilities  and  Advantages 

In  view  of  the  limitations  encountered  in  attempting  to  evaluate  the  pulse* 
pair  estimators  by  means  of  real  Porcupine  Radar  weather  return  data*  (see 
3.  1  above)  more  extensive  experimental  verification  of  their  theoretically 
predicted  performance  was  undertaken  on  the  present  contract  using  simulated 
radar  weather  return. 

The  use  of  simulation  techniques  to  generate  the  input  data  for  evaluating 
the  performance  of  the  estimators  provides  many  advantages  which  overcome  the 
principal  limitations  of  the  use  of  real  radar  data  for  evaluation  purposes.  Pri¬ 
marily,  it  makes  possible  the  generation  of  statistically  significant  numbers  of  T 

simulated  complex  radar  samples  with  prescribed  spectral  distributions  (mean 
frequency  and  width)and  S/N.  This  permits  the  generation  of  data  tapes  with 
controllable  parameters  made  to  order  and  in  sufficient  quantities  per  case  to 
suit  the  evaluation  needs.  Since  the  true  parameters  are  known  in  advance,  a 
direct  and  precise  determination  of  estimator  performance  accuracy  can  be  made 
at  all  times.  The  ensemble  of  precise  results  can  then  be  analyzed  by  the  com¬ 
puter  to  determine  the  statistical  behavior  of  the  estimators  against  the  synthetic 
distributed  weather  targets  as  a  function  of  the  controllable  parameters,  number 
of  samples  entering  into  the  estimate,  etc. ,  and  the  results  compared  with  the 
theory. 

The  particular  cases  chosen  for  the  simulation  tests  were  tailored  to  facilitate 
comparison  with  the  theoretical  performance  predictions  derived  in  Section  2. 

A  PRF  of  3300,  the  highest  available  for  the  Porcupine  Radar  was  used  for  all 
cases.  A  spectral  mean  of  300  Hz  and  width  of  78  Hz  were  selected  for  the 
initial  cases.  These  values  are  equivalent  to  a  weather.target  velocity  mean 
and  spread  of  approximately  8  meters  per  second  and  2  meters  per  second 
respectively,  as  observed  with  the  Porcupine  Radar.  To  vary  the  S/N  para¬ 
meter,  two  extreme  values  of  S/N  were  chosen,  including  0  dB  and  oodB  (no 
noise  present)  and  two  cases  closer  to  the  poor  S/N  condition.  To  conveniently 
vary  the  theoretically  significant  wT  parameter,  several  multiples  of  the  initial 
mean  and  spread  were  selected  to  yield  target  velocity  distributions  with  16  meters 
mean,  4  meters  spread;  and  24  meters  mean,  6  meters  spread.  Each  of  these 
cases  was  varied  over  the  range  of  four  S/N  values  established  initially. 
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An  overall  block  diagram  of  the  software  simulation  system  utilised  for 
investigating  the  pulse  pair  technique  appears  in  Figure  9*  The  system 
is  capable  of  generating  sequences  of  simulated  radar  doppler  return  signals 
with  controllable  parameters;  estimating  the  mean  and  width  of  the  spectral  dis¬ 
tributions  by  both  pulse  pair  and  spectral  analysis  teenniques;  performing  statis¬ 
tical  analyses  of  estimator  accuracy;  and  plotting  and/or  printing  out  the  results. 

The  simulation  system  is  programmed  in  FORTRAN  and  was  run  on  the 
CDC  6600  Computer  at  AFCRL«.  To  conveniently  suit  the  capacity  and  capabilities 
of  the  computer,  the  simulation  system  was  divided  into  separable  steps,  which 
could  be  run  at  different  times.  To  facilitate  such  division  into  steps,  the  se¬ 
quences  of  simulated  radar  doppler  return  which  were  to  be  used  as  input  data  to 
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the  signal  processors  under  investigation  were  recorded  on  magnetic  tape  for 
later  processing  and  reference.  A  library  of  such  input  data  tapes  was  generated 
to  provide  desired  numbers  of  statistically  independent  sequences  of  simulated 
return,  with  desired  combinations  of  spectral  parameters  and  S/N.  The  signal 
processing  schemes,  and  modifications  thereto,  could  thereby  be  statistically 
evaluated  and  compared  on  the  basis  of  identical  input  data.  Provision  was  also 
made  to  visually  verify  the  validity  of  the  simulated  data;  designated  sequences 
of  returns  from  among  those  generated  could  be  selected  for  spectral  analysis 
by  conventional  FFT  and  automatically  plotted  for  visual  check. 

The  forthcoming  sub-sections  deal  in  greater  detail  with  the  generation  of 
the  simulated  doppler  returns;  the  comparative  performance  of  the  pulse  pair 
and  spectrum  analysis  estimators  against  the  simulated  input  data;  and  a  com¬ 
parison  of  pulse  pair  estimator  performance  as  experimentally  determined  by 
simulation  with  that  predicted  f'om  theory.  Excellent  correlation  between  theory 
and  experiment  are  shown. 

3.2.2  Generating  Sequences  of  Simulated  Doppler  Radar  Target  Return 
with  Controllable  Parameters 

In  the  real  Porcupine  Radar  case,  a  uniform  PRF  was  uBed  to  generate 
sequences  of  uniformly  spaced  coherent  echoes  reflected  from  distributed 
meteorological  targets  as  observed  in  a  fixed  range  gate.  The  uniformly  spaced 
echoes,  in  the  form  of  digitized  complex  numbers  derived  from  the  phase  detected 
coherent  video  output, were  then  processed  two  at  a  time  to  yield  the  pulse  pair 
estimates  of  the  target  spectral  parameters. 

The  object  of  the  simulation  program  is  to  generate  similar  sequences  of 
uniformally  spaced  complex  numbers  equivalent  to  radar  echoes  reflected  from  a 
distributed  target  with  a  given  radial  velocity  distribution,  as  sampled  at  a  given 
uniform  PRF  and  with  a  given  output  S/N.  The  problem  is  mathematically  equiva¬ 
lent  to  that  of  numerically  generating  a  sequence  of  sampled  values  of  a  stationary 
Gaussian  process  having  a  specified  correlation  function  or  power  spectrum.  The 


method  employed  is  adapted  from  Levin 
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A.  Theory 

The  following  parameters  of  the  desired  sequence  must  be  controllable: 


Term 

Symbol 

Unit: 

MEAN  DOPPLER 

f 

o 

Hz 

SECOND  MOMENT 

2 

m2  =  w 

Hz* 

SAMPLING  RATE 

PRF 

Hz 

S/N  RATIO  (per  pulse) 

B 

dB 

1.  Spectral  Distribution  -  An  approximation  to  the  desired  spectral 
distribution  of  the  radar  video  data  is  achieved  by  playing  Gaussian 
noise  of  zero  mean  and  unity  variance  through  a  filter  consisting  of 
two  cascaded,  low -pass  networks.  The  transfer  function  of  the 
filter  is  given  by: 


H(s)  = - -  where  a  =  2nw  and  (3-5) 

(a  +  a)2  a  =  complex  frequency  (rad/sec) 

From  this  expression  it  is.possible  to  derive  the  second  moment  and 
half  power  point  of  the  filter.  The  second  moment  is  given  by: 

2  +j® 

H(s)  H(-s)(-si)  ds  (3-6) 


“2 


M 


=  w 


H(s)  H(-s)  ds 


The  half  power  point  (-3dB)  occurs  at 


W3=  a  [  nT  2-1  ] 


Vz 


(3-7) 
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By  ^-transform  theory,  the  recursion  relation  for  the  equivalent  digital 
filter  is  given  by  the  impulse  invariant  transformation;  namely,  if 


H(s)  = 


(fl+a) 


then, 


„-i _  2  ~ 

Z  a  1  e 

H(Z)  =  “  "  •  where  T  * 

a  .  e*wTz-‘  )2  PRF 


(3-8) 


The  resulting  recursion  relation  is  then 

z  ^  -aT,  ,  -  aT  -2aT 

y  =  x  y  T  e  +  2  e  y  -  e  y 
n  n-l  n-l  n-c 


(3-9) 


Simulation  of  radar  quadrature  data  is  achieved  as  illustrated  in 

Figure  10.  The  sequences  fxj.  are  two  independent  Gaussian  sequences 
which  are  then  played  through  the  same  filter,  H(Z),  specified  by  (3-8)  and 


(3-9). 

Complex  multiplication  of  the  output  sequences  (y,.  y2]  by  e 
shifts  the  mean  frequency  fromzero  to  f  . 


j2uf  ^nT 


Transient  effects  of  the  filter  are  minimized  by  discarding  a  nurrber 
of  initial  terms.  This  number  has  been  set  at  three  times  the  fractional 
bandwidth  of  the  filter  (with  respect  to  the  PRF)  and  it  is  given  by; 


3  •  (PRF) v 

No.  of  discarded  initial  terms  =  - 

w 


(3-10) 


2.  Signal  to  Noise  Ratio  -  When  Gaussian  noise  is  added  to  sequences 

(y1  ,  y'  ]  (See  Figure  10  ),  the  per-pulse  signal  to  noise  ratio  (S/N)  is 
i  i 

specified  by: 
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S/N  = 


a4  T*  (3e 


-aT 


+  e 


‘n(e  aT-  * 


ns — '  ■- B  tdB) 


where  a^is  the  variance  of  the  added,  zero  mean,  Gaussian  noise.  A  given 

is 

-B/10 


S/N  ratio  is  then  achieved  by  specifying  a^as 


*  -  (3e'aT+  eaT)  10 

°N =  aT  -  aT  . 


(3-12) 


(z 


-  e 


where  B  is  the  desired  S/N  ratio  in  dB. 


Mathematical  derivations  relating  to  the  above  processes  appear  in 
Appendix  B. 

B.  Description  of  Computer  Program 

The  program  developed  for  generating  the  simulated  doppler  target  return 
on  the  AFCRL.  CDC  6600  is  called  DOPGEN.  The  flowchart  for  DOPGEN  is  shown 
in  Figure  11  .  and  the  program  statements  coded  in  FORTRAN  appear  in 
Appendix  Cl. 

In*put  to  the  program  is  in  the  form  of  data  cards  on  which  are  specified: 

(1)  mean  Doppler  frequency; 

(2)  standard  deviation  of  Doppler  frequency; 

(3)  per-pulse  signal-to-noise  ratio; 

(4)  FFT  plot  control;  and, 

(5)  number  of  blocks  (N)  of  1024  complex  samples. 

Upon  reading  the  first  data  card,  the  program  generates  1024  complex  samples 
which  simulate  Doppler  return  (a  PRF  of  3300  is  used  throughout).  The 
data  block  is  then  written  on  magnetic  tape  along  with  twelve  ancillary  words 
containing  information  about  the  data.  This  procedure  continues  until  N  inde¬ 
pendent  data  blocks  are  written  with  the  specified  Doppler  mean  and  standard 
deviation.*  A  second  data  card  can  be  used  to  change  the  mean  and  the  width 
of  the  Doppler  spectrum  or  to  change  the  S/N  ratio.  If  no  more  data  blocks 

*  J>ince  each  data  block  occupies  about  3.5  ft.  of  magnetic  tape,  685  data  blocks  will 
fit  on  a  standard  2400  ft.  long  tape. 


3-  14 


Complex  Samples 


Figure  11.  Flowchart  of  DOPGEN 


are  to  be  generated,  as  indicated  by  a  blank  data  card,  the  program  rewinds 
the  tape  and  reads  each  block  to  verify  the  readability  of  the  tape.  In  addition, 
for  those  data  blocks  on  which  the  FFT  plot  control  parameter  (one  of  the 
ancillary  words)  is  set  to  one,  the  program  performs  an  FFT  on  the  data  and 
plots  the  transformed  data  on  a  Calcomp  Plotter  (available  at  AFCRL).  This 
is  used  as  a  check  to  verify  the  validity  of  the  simulated  data  -  a  sample  plot 
is  shown  in  Figure  12. 

A  library  of  390  separate  sequences  of  1024  complex  samples  was  generated 
by  the  DOPGEN  simulation  program  and  recorded  on  magnetic  tape  for  use  au  input 
data.  Each  sequence  of  1024  samples  constitutes  a  data  block.  Three  pairs  of 
mean  doppler  and  doppler  standard  deviation  were  run,  each  at  four  conditions 
of  S/N;  a  PRF  of  3300  pps,  representative  of  the  Porcupine  Radar,  was  used  in 
all  cases.  Thirty  independent  data  blocks  with  specific  parameters  as  shown  in 
Table  3  were  generated  for  each  combination  of  parameters.  A  sample 
spectral  analysis  plot  for  each  of  the  twelve  combinations  is  presented  in 
Appendix  E. 


Block  No. 

Mean 

Doppler 

(Hz) 

Std.  Dev. 
Doppler 
(Hz) 

S/N 

(dB) 

Tape 

41-70 

300 

78 

5 

RPK07 

101-130 

600 

156 

00 

RPK08 

131-160 

600 

156 

0 

RPK08 

161-190 

600 

156 

5 

RPK08 

191-220 

600 

156 

15 

RPK08 

221-250 

900 

234 

00 

RPK08 

251-280 

900 

234 

0 

RPK08 

251-280 

900 

234 

5 

RPK08 

311-340 

900 

234 

15 

RPK08 

341-370 

300 

78 

oo 

RPK09 

371-400 

300 

78 

0 

RPK09 

401-430 

300 

78 

5 

RPK09 

431-460 

300 

78 

15 

RPK09 

Table  3.  List  of  Recorded  Sequences  of  Simulated  Doppler  Return 
Generated,  with  Parameters  as  shown. 

(PRF  =  3300  pps  in  all  cases) 
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3.2.3  Performance  of  Pulse  Pair  Estimators 


A.  Computer  Program 

As  shown  In  Figure  9  of  the  overall  simulation  system,  the  tape 
recorded  sequences  of  simulated  doppler  return  are  processed  by  the  pulse 
pair  algorithms  to  yield  estimates  of  spectral  mean  and  width  as  a  function  of 
the  number  of  pulse  pairs  entering  into  the  estimate.  The  resulting  estimates 
are  then  compared  with  the  known  given  spectral  parameters  of  the  input  data 
and  the  performance  of  the  pulse  pair  estimators  is  statistically  analyzed  and 
summarized  to  produce  printouts  of  estimator  accuracy  versus  number  of  pulse 
pair  samples  for  each  simulated  test  condition. 

The  program  for  pulse  pair  processing  and  statistical  evaluation  is  called 
PPSTAT.  A  flowchart  of  PPSTAT  appears  in  Figure  13,  and  the  coded 
FORTRAN  statements  for  this  program  may  be  found  in  Appendix  C-2.  As 
shown  in  the  flowchart,  before  reading  the  magnetic  tape  generated  by  DOPGEN, 
PPSTAT  reads  from  a  data  card  two  data  block  numbers  (the  numbers  are  assigned 
by  DOPGEN  and  are  contained  in  the  ancillary  data  associated  with  each  block). 
Pulse  Pair  estimates  are  then  made  on  consecutive  data  blocks  starting  with  the 
first  block  number  and  continuing  to  the  data  block  with  the  second  block  number. 
The  estimates  are  made  on  samples  of  size  16,  32,  64,  128,  256,  512,  and  1024. 
After  operating  on  the  variable  sample  sizes  from  each  data  block,  the  program 
performs  a  statistical  summary  of  all  the  Pulse  Pair  estimates  to  yield  estimator 
accuracy  versus  sample  size. 

B.  Results 

The  hundreds  of  pulse  pair  estimates  of  spectral  mean  and  width  produced, 
analyzed,  and  printed  out  by  the  simulation  system  are  presented  in  a  manner 
intended  to  facilitate  comparison  with  theoretically  predicted  performance  capa¬ 
bilities.  The  close  correlation  shown  between  theory  and  experiment  in  the 
1  Figures  that  follow  confirms  the  validity  of  the  theory  under  the  practical  condi¬ 
tions  of  noise  background,  correlation  between  pulse  pairs,  etc.  Analytic  predic¬ 
tions  of  pulse  pair  estimator  performance  for  particular  sets  sf  conditions  can, 
therefore,  be  made  with  confidence. 

The  following  comparisons  between  pulse  pair  experimental  data  and  theory 
have  been  calculated  and  plotted  from  the  results  of  the  simulation  tests. 
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Figure  13,  Flowchart  of  PPSTAT 
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Accuracy  of  Mean  Frequency  Estimate*  vs  Number  of  Pulse  Pairs 

Accuracy  of  Mean  Frequency  Estimates  vs  Spectral  Width 

Accuracy  of  Spectral  Width  Estimates  vs  Number  of  Pulse  Pairs 

Accuracy  of  Spectral  Width  Estimates  vs  Actual  Spectral  Width 

Accuracy  of  Spectral  Width  Estimates  vs  Signal-to»Noise  Ratio 
(for  fc  s  300  Hi;w  a  7  8  Hs) 

Accuracy  of  Spectral  Width  Estimates  vs  Signal-to-Noise  Ratio 
( for  f  =  600  Hz;  w  =  156  Hz) 

Accuracy  of  Spectral  Width  Estimates  vs  Signal-to-Noise  Ratio 
(for  f c  =  900  Hz;  w  =  234  Hz) 
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Figure  14.  Theoretical  Accuracy  of  Pulse  Pair  Estimates  of 
Mean  Frequency  (f0)  vs  Number  of  Pulse  Pairs 
Compared  with  Similation  Results 
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Figure  20.  Theoretical  Accuracy  of  I'uhe  Pair  Estimates  of  Spectral 
Width  (w)  vs  S/N  Compared  with  Simulation  Results  (for 
f  =  900  Hz;  w  =  234  Hz) 


3.2.4  Performance  of  Spectral  Analysis  Estimator 


A.  Computer  Program 


In  order  to  provide  a  standard  against  which  to  compare  the  performance  of 
the  pulse  pair  estimators  under  varying  conditions  of  signal  to  noise,  the  simulated 
doppler  return  generated  by  DOPGEN  was  also  analyzed  by  the  conventional  spectral 
analysis  method.  As  in  the  case  of  real  Porcupine  radar  data  (Section  3.  1.  1  above), 
the  spectral  parameter  estimates  were  obtained  in  two  steps.  First,  a  spectral 
analysis  is  performed  on  an  entire  data  block  to  yield  its  spectrum;  then  the  mean 
and  spectral  width  parameters  are  calculated  from  all  spectral  components  above  a 
given  threshold  level,  for  a  range  of  threshold  levels.  The  resulting  estimates  are 
printed  out  as  functions  of  threshold  level.  This  spectral  analysis  program  is 
designated  SPEC.  As  shown  in  the  flowchart.  Figure  21,  the  Spectral  Analysis 
program  searches  the  tape  generated  by  DOPGEN  and  reads  the  data  block  whose 
number  is  specified  on  the  input  data  card.  The  1024  complex  samples  are  then 
weighted  using  60  dB  Dolph-Chebyshev  weighting  coefficients  and  input  into  the  FFT 
sub-routine.  If  desired,  the  FFT  output  will  be  plotted.  Next,  the  mean  aud 
standard  deviation  are  calculated  ba6ed  on  all  FFT  output  points  above  a  specified 
threshold  level.  The  threshold  level  is  automatically  incremented  in  steps  of  2  dB 
starting  at  2  dB  below  the  maximum  FFT  output  point  and  continuing  to  40  dB  below 
the  maximum  point.  The  mean  and  standard  deviation  values  are  printed  on  the 
line  printer  as  a  function  of  threshold  level.  The  above  procedure  is  repeated  for 
each  input  data  card  until  a  blank  card  is  encountered.  SPEC  coded  in  FORTRAN 
appears  in  Appendix  C3. 

B.  Results 

It  was  concluded  in  the  discussion  of  results  with  Porcupine  Radar  input  data 

A  t 

(Section  3.  1  above)  that  the  spectral  analysis  estimator  of  fQ  was  definitely  satis¬ 
factory.  Therefore,  no  further  analysis  and  evaluation  of  the  performance  of  this 
estimator  was  done  with  the  results  of  the  simulation  tests. 

On  the  other  hand,  there  was  interest  in  how  well  the  threshoided  spectral 
analysis  estimator  of  spectral  width  (w)  would  show  up  against  the  simulated  input 
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Figure  21.  Flowchart  of  SPEC 
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Typical  results  are  plotted  in  Figures  22,  23,  and  24.  Each  figure 
depicts  the  width  estimates  produced  by  a  single  block  of  samples  for 
different  combinationo  of  input  data  parameters,  fQ,  w  and  S/N  ,  as  a 
function  of  threshold  level.  As  was  indicated  earlier,  no  claim  is  made 
that  the  thresholdec  estimator  of  Eq.(3-4)  is  optimum.  It  is,  however, 
straightforward  and  readily  implementable  by  computer.  To  preclude 
possible  bias  effects  from  aliasing,  the  spectral  width  was  calculated  about 
the  true  mean,  fQ.  Despite  these  precautions,  the  width  estimates  again  are 
shown  to  be  very  sensitive  to  the  value  of  the  threshold.  Only  in  the  cases 
of  very  high  S/N  do  the  estimates  asymptotically  approach  the  true  width. 


Figure  22.  Spectral  Analysis  Width  Estimates  (w)  vs 
Threshold  Level  for  Simulated  Input  Data 
(f  =  300  Hz;  w  =  78  Hz) 
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Figure  24.  Spectral  Analysis  Width  Estimates  (w)  vs 
Threshold  Level  for  Simulated  Input  Data 
(f  =  900  Hz:  w  =  234  Hz) 


3.2.5  Comparison  of  the  Width  Estimators 


The  accuracy  achieved  by  the  pulse  pair  width  estimator  in  the  simulation 
tests  is  compared  in  Table  4  with  the  accuracy  predicted  from  theory,  and  the 
estimates  yielded  by  the  spectral  analysis  estimator.  The  results  are  shown 
for  all  12  cases  of  fQ,  w  and  S/N.  (It  should  be  noted  that  the  pulse  pair  data  is 
presented  in  the  form  of  the  standard  deviation  of  the  estimates  from  the  true 
value,  based  on  a  statistical  analysis  of  20  estimates  on  as  many  independent 
blocks  of  simulated  data  for  each  case  shown;  whereas  the  spectrum  analysis 
results  are  for  a  single  block  of  samples  for  each  case,  as  a  function  of  threshold 
level. ) 


Table  4  and  Figures  18-20  show  that  the  pul6e  pair  width  estimates  obtained 
with  the  1024  samples  and  S/N  in  excess  of  5  dB  per  pulse  produces  width  errors 
of  less  than  3  dB  for  all  cases  tested.  Translated  into  meteorological  terms,  for 
the  Porcupine  radar,  actual  widths  of  .  2  m/ sec  were  estimated  at  .  3  m/ sec 
with  the  equivalent  at  a  1/3  sec  look  (at  3300  pulses/sec)  at  a  signal -to -noise  ratio 
of  5  dB. 

For  a  signal  which  is  a  Gaussian  random  variate  with  spectral  density  function 


!• 


frr™ vT 


imbedded  in  white  noise  with  spectral  density  coefficient  Nq,  the  ratio  of  the  peak 
spectral  value  to  the  rms  noise  background  level  in  the  frequency  domain  is  given  by 


r  .  SisL  ,  ’ 

'  N  7^^ 


If  this  datum  is  sampled  at  the  rate  1/T,  the  per  pulse  eignal-to-noise  ratio  is 
given  by 


c(f)  df 
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Table  4.  Comparison  Between  Pulse  Pair  &  Spectrum  Analysis;  Estimates  of  Spectral  Width 
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Thus  the  spectral  peak  output  of  the  simulated  signals  exceeds  the  input  sign&l-to* 
noise  ratio  by  approximately  the  values  shown  in  Table  5. 

Table  5 

Peak  Sisnal-to-Noise  Improvement  of  Spectrum  Analyzer 


r  _ 

W  (Hz) 

D  in  dB 

i  * 

78 

12.  3 

\ 

r 

156 

9.  3 

i 

234 

6.  3 

This  calculation  shows  that  the  spectrum  analyzer  output  reduction  to  mean 
and  variance  may  be  expected  to  perform  extremely  poorly  when  the  threshold  is 
set  at  approximately  D  +  S/N  (in  dB)  below  the  peak  output.  These  points  corres¬ 
pond  approximately  to  the  break  points  shown  on  Figures  22  to  24.  Assuming  that 
the  spectrum  analyzer  threshold  had  been  set  5  dB  above  the  noise  background  level 
(high  enough  to  assure  that  the  noise  peaks  would  not  disturb  the  measurement). 
Figures  22  to  24  indicate  that  the  measured  widths  would  be  as  shown  in  Table  6. 

It  should  again  be  noted  that  the  data  shown  on  Figures  22  -  24  represents 
only  a  single  spectral  sample  and  that  the  thresholds  shown  are  measured  with 
respect  to  the  sample  peak  value.  Table  6  attempts  to  incorporate  this  factor  by 
using  the  breakpoints  of  the  measurement,  as  the  threshold  level.  These  are 
at  approximately  the  threshold  shown  in  Table  6. 

In  summary,  the  spectrum  analyzer  measurements  appear  to  consistently 
underestimate  the  width  parameter  while  the  pulse  pair  technique  overestimates 
it.  The  simulation  indicates  that  the  accuracies  available  for  basically  the  same 
data  conditions  are  roughly  comparable  for  the  cases  tested.  Figures  22  -  24 
indicate  a  marked  sensitivity  to  threshold  settings  which  is  not  directly  evident  in 
the  pulse  pair  measurements.  It  is  important  to  note  that  the  pulse  pair  measure¬ 
ment  of  width  uses  an  estimate  of  the  noise  background  level  also  and  thus  has  a 
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Comparison  of  Pule*  Pair  and 
Spectrum  Analyser  Derived  Width  Estimates 

(Spectrum  Analyser  Threshold  Set  5  dB  Above  Average  Background  Neve  el 
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corresponding  parameter  uncertainty  in  its  process.  In  other  words,  both 
techniques  require  a  measurement  of  the  noise  alone  background  to  obtain 
useful  final  results. 

The  conclusion  to  be  drawn  from  this  comparison  is  that  the  pulse  pair 
technique  shows  no  evident  lack  of  sensitivity  as  compared  with  the  spectrum 
analyzer  approach  in  situations  characteristic  of  the  meteorological  applications. 
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4.  0  HARDWARE  REALIZATION 


This  section  describee  a  real  time  digital  processor  design  to  reduce  the 
pulse  pair  technique  to  practice.  A  brief  theoretical  review  of  the  signal  processing 
algorithms  to  be  implemented,  is  followed  by  discussions  and  block  diagrams  of  the 
processor  and  its  key  components.  Design  goals  for  the  proposed  pulse  pair  pro¬ 
cessor  are  defined,  and  its  performance  compared  with  CMF,  Analog  Spectrum 
Analyzer,  and  Fast  Fourier  Transform  processors.  Finally,  the  design  tradeoffs 
to  be  considered  during  the  design  phase  are  discussed. 

4.  1  General  Description  of  the  Pulse-Pair  Processor 

This  processor  operates  on  the  output  of  a  radar,  and  is  capable  of  pro¬ 
ducing  real-time  independent  estimates  of  the  spectral  mean  and  width  of  a 
Doppler  radar  signal  for  each  range  cell  viewed  by  the  radar. 

The  processor  accepts  the  radar  IF  signal  and  system  trigger  as  inputs 
and  computes  the  mean  frequency  (mean  velocity) -and  spectral  width  (velocity 
distribution)  in  each  of  256,  512,  768  or  1024  range  cells,  selectable  by  a  front 
panel  switch.  Each  range  cell  has  a  width  of  0.  5,  1  or  2  microseconds  selectable 
by  a  front  panel  switch.  The  processor  produces  both  analog  outputs  for  viewing 
on  real  time  displays,  and  digital  outputs  which  may  be  remoted  or  recorded 
for  future  computer  processing.  An  internal  digital  automatic  gain  control  sys¬ 
tem  independently  measures  the  signal  amplitude  and  adjusts  the  gain  of  the 
processor  for  each  range  cell  processed  to  maintain  a  90  dB  dynamic  range. 

Although  the  pulse  pair  technique  will  work  with  either  coherent  or  non¬ 
coherent  radars,  a  coherent  receiver  output  is  necessary  in  order  to  provide 
complex  samples  of  the  target  return  to  the  processor.  In  addition,  to  achieve 
independent  spectral  estimates  from  range  cell  to  range  cell,  the  receiver  gain 
must  have  a  wide  dynamic  range  and  be  range  cell-to-cell  controllable.  Accord¬ 
ingly,  the  receiver  requirements  must  be  included  as  part  of  the  processor  con¬ 
siderations;  how  much  of  an  existing  radar  receiver  can  be  used  as  is  and  bow 
much  must  be  modified  or  incorporated  as  part  of  the  processor  depends  upon 
the  characteristics  of  the  particular  radar  whose  output  is  to  be  processed.  The 
receiver  for  the  Pulse  Pair  Processor  described  below  is  interfaced  to  match  the 
Porcupine  Radar. 
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4.  2  Technical  Description  of  the  Pulse  Pair  Processor 


The  PPP  performs  spectral  mean  and  width  estimations  according  to  equa¬ 
tions  (4-1)  and  {4-3)  respectively,  which  are  equivalent  to  earlier  equations  (2-9) 
and  (2-10).  They  are  rewritten  here  in  more  convenient  form  for  the  discussion 
of  hardware. 
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where  the  z's  are  complex  samples  of  signal  plc_  noise,  and  Nq  is  the  receiver 
noise  power  in  the  absence  of  signal,  z*  denotes  the  complex  conjugate  of  z. 
Equation  (4-  1)  estimates  mean  frequency  by  determining  the  phase  differences 
in  pairs  of  complex  samples  (zk+1»  zk>»  averaging  these  differences  to  find 
the  mean  phase  difference,  and  scaling  this  mean  phase  difference  by  the  PRF 
to  form  a  mean  frequency  estimate.  The  phase  differences  averaged  are 
weighted  according  to  sample  amplitude  so  that  phase  differences  from  strong 
returns  count  more  than  differences  from  weak  returns.  The  results  in  Section 
3  showed  that  as  the  number  of  samples  processed  increases,  the  accuracy  of 
fQ  improves  asymptotically. 

A 

Equation  (4-2)  estimates  spectral  width  w  in  terms  of  signal  coherence, 
including  a  correction  for  noise,  Nq, 
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which  may  be  measured.  The  parameter  a  is  well-defined  for  all  sampled 
signals  including  those  with  multimodal  spectra.  It  may  be  seen  from  equation 
(4-3)  that  in  the  absence  of  noise  a  =  1  for  sinusoidal  signals,  and  that  a  does  not 
depend  on  mean  signal  frequency,  which  may  exceed  the  complex  sampling  fre¬ 
quency  (PRF).  As  is  usual  in  sampled  data  systems,  the  signal  bandwidth  must 
not  exceed  the  PRF. 
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Figure  25  shows  the  signal  flow  of  the  PPP  algorithm  applied  to  the 
sampled  Doppler  signal  from  one  range  cell.  A  delay  and  multiplier  are  used 
to  form  the  products  z*  ,  which  are  then  averaged.  The  sample  energies^ 

{  2  |  2  are  also  averaged,  and  function  generators  are  used  to  compute  fQ  and  w. 

There  are  a  number  of  alternate  ways  of  deriving  the  noise  correction  term 
(N2  )  required  for  the  calculation  of  the  width  estimate.  One  such  method  is  de¬ 
picted  in  Figure  25.  A  measurement  of  receiver  noise  (nQ)  (which  may  be  made, 
for  example,  at  the  output  of  the  receiver  at  a  range  in  excess  of  that  in  which  any 
signal  return  may  be  expected)  is  corrected  for  the  receiver  gain  applicable  to  the 
particular  range  cell,  and  introduced  as  the  noise  correction  (Nq)  in  the  width  esti 
mate  (w).  The  correction  factor  is  obtained  by  calibrating  the  receiver.  The  mean  j 

power  output  (p),  which  is  made  up  of  averaged  measurements  of  signals  plus  noise, 
is  delivered  for  use  as  the  range  cell's  control  signal  to  set  the  receiver  ACC  at  j 

its  range.  t 

j 

Detailed  design  of  the  noise  correction  scheme  has  been  deferred  pending 

i 

a  hardware  tradeoff  study  of  alternate  schemes,  to  be  performed  as  part  of  a  \ 

future  hardware  fabrication  effort.  The  description  of  the  system  design  which 
follows  does  not  include  the  noise  correction  feature. 

4.  3  System  Design 

The  PPP  system,  shown  in  Figure  26,  comprises  a  fast  AGC  IF  amplifier, 
a  two-point  correlator,  three  (or  four)  digital  integrators,  a  timing  generator  and 
several  function  generators.  The  PPP  accepts  30  MHz  IF  input  and  produces 

A  A 

digital  and  analog  f  ,  w,  range  and  integrated  video  outputs.  The  major  sub¬ 
systems  of  the  PPP  are  scribed  below. 

4.3.1  IF  Amplifier  Phase  Detector 

Achievement  of  the  processor  90  dB  dynamic  range  requires  independent 
gain  control  in  each  range  cell  processed,  which  implies  that  the  AGC  voltage 
must  be  a  time  varying  signal  with  a  bandwidth  comparable  to  the  IF  bandwidth. 

Such  a  system  requires  an  IF  amplifier  of  unusual  design  and  closely  integrated 
IF  amplifier  and  AGC  loop  designs.  For  this  reason,  the  IF  amplifier  and  phase 
detectors  are  included  in  the  PPP. 
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An  IF  amplifier-phase  detector.  Figure  27,  was  recently  designed  for  the 
PORCUPINE  radar  which  is  acceptable  for  the  PPP  with  slight  modifications.  The 
amplifier  achieved  approximately  60  dB  ACC  control  range,  with  AGC  applied  to 
two  out  of  a  total  of  three  identical  gain  stages.  ACC  response  time  was  less  than 
200  ns.  By  adding  AGC  to  the  third  gain  stage,  this  amplifier  is  expected  to 
achieve  a  90  dB  AGC  control  range  without  degradation  of  response  time. 

4.  3.  2  Dipitizer 

The  digitizer  accepts  in-phase  and  quadrature  bipolar  video  and  produces 
digitized  estimates  of  these  inputs.  Eight  bit  quantization  appears  more  than 
adequate  to  compute  f  to  the  accuracy  of  the  algorithm.  A  single  A/D  converter, 
such  as  the  ILC  Data  Devices  VADC-B,  may  be  switched  to  quantize  both  input 
channels.  This  is  preferable  to  the  use  of  two  simultaneous  A/D  converters 
because  the  single  converter  is  less  expensive  and  because  problems  of  balancing 
the  input  channels  are  simplified  by  the  use  of  a  single  A/D. 

The  suggested  A/D  converter  can  take  an  analog  sample  every  150  ns  and  is 
compatible  with  a  range  resolution  of  500  ns.  A  150  ns  delay  line  must  be  used  in 
one  input  channel  to  compensate  for  the  A/D  sampling  delay.  This  A/D  includes 
a  sample  and  hold  circuit,  operates  from  standard  +5V  and  ±  15V  power  supplies 
and  is  contained  on  a  single  8.8"  x  4.  5"  circuit  card. 

4.  3.  3  Correlator 

The  correlator  accepts  complex  video  samples  z^  from  the  digitizer  and 
generates  products  of  the  form  z^  ^  zj£  where  the  sample  z^j  i®  delayed  by  one 
range  sweep  time  from  the  sample  z^.  This  delay  is  accomplished  using  MOS 
dynamic  shift  registers  organized  as  16  bits  x  512  words.  Recently  available 
1024  bit  registers,  6uch  as  the  Intel  1404A,  allow  this  delay  to  be  fabricated 
quite  compactly. 

The  complex  multiplier  which  forms  the  products  will  be  similar 

to  the  multiplier  developed  for  a  proprietary  Raytheon  FFT  pulse  compression 
system.  This  multiplier,  shown  in  Figure  28,  employs  programmed  arrays  and 
adders  to  perform  a  24  x  24  bit  complex  multiply  and  add  every  0.  5  P-s.  Approxi¬ 
mately  280  integrated  circuits  (DIP'S)  are  required,  hence,  the  multiplier  could 
be  fabricated  on  four  72-DIP  Augat  panels.  The  PPP  requires  only  a  16  x  16  bit 
complex  multiplier;  one  circuit  card  could  probably  be  saved  by  redesigning  the 
FFT  multiplier  for  the  less  demanding  PPP  requirement. 
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4.  3.  4  Digital  Integrators 


A  minimum  of  three  digital  integrators  are  required,  two  to  compute  the 
real  and  imaginary  parts  of  u  and  one  to  compute  average  power  p.  These  three 
integrators  will  be  similar  to  an  integrator  recently  designed  for  AFCRL,  to 
minimize  design  costs,  but  the  number  of  sweeps  integrated  will  only  be  variable 
up  to  1024  by  factors  of  2. 


An  integrated  power  measurement  is  required  to  operate  the  AGC,  This 
output  may  be  derived  either  from  the  power  channel  (p )  of  the  PPP  or  through  a 
separate  integrator.  The  integrator  which  produces  AGC  output,  whether  separate 
or  power  channel,  must  be  of  the  exponential  rather  than  square  moving  window 
type.  Theory  shows  that  the  gain  and  phase  of  the  square  window  integrator. 


/»■■■■ 

I  gy^l  3  ,J3>F~  VI-  cos  tuT 


(4-4) 


4»(<u)  =  Tan'1^'  1 


sin  u)J 
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makes  it  unsuitable  for  use  in  stable  feedback  loops  with  loop  gain  greater  than  1. 

A  design  study  should  determine  the  performance  differences  in  PPP' a  employing 
only  three  exponential  window  integrators,and  PPP'e  employing  three  square 
window  integrators  to  compute  u  and  p  and  one  exponential  integrator  for  AGC, 

4.  3.  5  Function  Generators 

Three  function  generators  are  required  for  the  PPP:  a  squared  magnitude 
generator,  a  log  generator  and  an  output  function  generator.  As  shown  in  Figure  27, 
all  function  generators  will  consist  of  programmed  arrays  combined  with  registers, 
gating  circuits  and  adders  appropriate  to  the  function  computed. 

Programmed  arrays  will  contain  the  functions  x2.  In  x,  tan“1e“X  and 

"X  l/Z  * 

(1  “  e  )  •  Each  programmed  array  will  consist  of  two  integrated  circuits 

®^ch  as  the  Monolithic  Memories  MM6300,  The  Mrv.^300  is  organized  as  an  8-bitin, 

4-bit  out  array  which  may  be  programmed  to  contain  any  desired  bit  pattern;  two 

such  circuits  form  an  8-bit  in,  8-hit  out  digital  function  generator.  Approximately 

50  ns  after  an  8-bit  input  word  is  supplied  to  the  two  programmed  arrays,  an  8-bit 

output  word  is  available,  allowing  each  pair  of  arrays  to  compute,  and  an  output 

register  to  store,  up  to  4  outputs  during  an  0.  5  us  range  cell. 

These  arrays  may  be  ordered  preprogrammed  by  the  manufacturer  or  may  be 
programmed  by  the  user.  Some  vendors  (e.g.  ,  Motorola)  sell  programmers  for 
their  arrays.  Raytheon  is  currently  constructing  a  programmer  for  the  Monolithic 
Memories  array  in  connection  with  a  high  speed  pulse  compression  FFT  develop¬ 
ment;  this  programmer  is  expected  to  be  available  for  use  in  PPP  fabrication. 
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4. 4  Performance 

The  detailed  performance  of  the  PPP  will  depend  upon  a  number  of  design 
decisions,  such  as  the  AGC  loop  characteristics  and  the  exact  form  of  the  function 
generators;  operational  considerations  such  as  the  signal  characteristics  and  the 
number  of  pulses  integrated;  and  the  basic  accuracy  of  the  pulse  pair  algorithm. 
Analysis  of  these  factors  should  be  made  in  a  design  study  but  is  beyond  the  scope 
of  this  hardware  feasibility  report. 

By  assuming  a  set  of  design  goals  for  PPP  hardware,  subject  to  change  in 
the  design  study,  a  few  general  comparisons  may  be  made  between  the  PPP  and 
other  related  signal  processing  hardware.  The  design  goals  in  Table 7  are  there¬ 
fore  assumed  for  the  purposes  of  comparison. 

4.4.1  Comparison  with  the  CMF 

The  Porcupine  Coherent  Memory  Filter  is  a  spectrum  analyzer,  and  does 
not  produce  estimates  of  spectral  mean  and  width  directly.  Modification  cf  the 
CMF  to  accommodate  different  pulse  repetition  rates  and  pulse  widths  is  expensive 
and  time  consuming.  By  contrast,  the  PPP  directly  produces  both  spectral  mean 
and  width  estimates  in  digital  form  in  real  time  and  operates  over  a  variety  of 
PRR'a  and  pulse  widths.  Its  quantitative  estimates  can  also  be  recorded  or 
remoted  for  further  processing,  plotting  or  display. 

Recently,  the  CMF  has  been  used  to  generate  displays  from  which  an  ob¬ 
server  can  estimate  mean  frequency.  The  PPP  can  generate  such  displays  with 
much  finer  time  and  frequency  resolution  by  addition  of  a  low  cost  voltage-to-delay 
converter.  In  addition,  the  PPP  mean  and  standard  deviation  outputs  can  be 
directly  displayed  on  a  PPI  scope  to  provide  Plan  Mean  Indicator  or  Plan  Width 
Indicator  displays.  These  displays  could  be  contoured  if  desired  by  using  existing 
contour  generator  equipment  (such  as  digital  integrators). 

The  following  compares  the  resolution  capabilities  of  the  CMF  and  PPP 


CMF 

PPP 

Unit 

Range  Resolution 

1860 

150 

Meters  (m) 

Range  Clements 

192 

1024 

- 

Velocity  Resolution 

0.  5  Xf  /50 
o 

0.  5  Xf  /  1024 
o' 

m/  s 

Dynamic  Range 

30 

90 

dB 

Input  Data  Rate 

1.75 

32 

M  Bits/ s 

Size 

1.0 

0.  3 

19"  rack  x  6' 

Table  7.  Proposed  Design  Goals  for  Pulse  Pair  Processor 


Signal  Input: 

Reference  Input: 

T rigger  Input: 

Analog  Outputs: 

Digital  Outputs: 

T rigger  Output: 

Range  Cells  Processed: 
PRF: 

Range  Cell  Width: 

Samples  Processed: 

Internal  Quantization: 

Dynamic  Range: 

AC.C: 

Input  Power: 

Mechanical  Configuration: 
Ambient  Temperature: 

C  onstruction: 


DESIGN  GOAL 

30  MHz  IF,  -80  dBm  to  +10  dBm 
Input  Impedance  50  ohms 

30  MHz,  1  v  p-p  into  100  ohms 

2  v  p  -  90  v  p  pulse 

0.  1  ps  -  100  ps  pulse  width 

Mean  frequency,  standard  deviation  mean  power 

and  bipolar  video 

Output  Impedance  -  50  ohms 

Mea:  .requency,  standard  deviation  mean  power 
and  cell  count  buffered 

TTL  level  outputs  10,  8,  and  10  bits  respectively 

Nominal  15  v  p  into  50  ohms  -  coincident  with 
1st  range  cell  output 

256,  512,  768  or  1024  nominal  selectable  by  front 
panel  switch 

300  pp's  -  5000  pps 

0.  5  ps,  1  ps,  or  2  pe  -  selectable  by  front  panel  switch 

16,  32,  64,  178,  256,  512,  or  1024,  selectable  by 
front  panel  switch 

Digital  signals  shall  be  represented  by  a  minimum 
of  8  bits  at  maximum  input  level 

Normal  and  quadrature  bipolar  video  shall  be 
represented  internally  by  a  minimum  o(  7  hits 
over  an  IF  input  range  of  100  pv  rms  -  1  v  rms 
with  1024  samples  processed 

The  processor  shall  include  a  dynamic  AGC  system 
to  provide  independent  gain  control  for  each  range 
ceil  processed. 

1  15  v  AC  +10C,  60  Hz  +2  Hz 

Mountable  in  a  19"  rack  cabinet 

o  o 

0  C  -  40  C 

According  co  best  commercial  practice 
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i  4. 4.  2  Comparison  with  Analog  Spectrum  Analyzers 

[  Presently  available  analog  spectrum  analyzers  such  as  the  Federal  Scientific 

|  "Ubiquitous"  may,  like  the  CMF,  be  used  to  generate  a  display  from  which  an 
/  observer  can  estimate  spectral  means  and  widths  indirectly.  Unlike  the  CMF  and 

|  PPP,  however,  such  analyzers  are  not  presently  capable  of  processing  the  entire 

S  output  of  a  weather  radar  in  real  time,  and  cannot  generate  any  real-time  displays. 
|  The  following  compares  the  resolution  of  the  PPP  and  Ubiquitous  analyzer. 


Ubiquitous 

PPP 

Unit 

Range  Resolution 

1200* 

150 

Meters  (m) 

Range  Elements 

10 

1024 

- 

Velocity  Resolution 

0.  S  Xf  /500 

0.  5  >  f  / 1024 
o 

m/  s 

Dynamic  Range 

40 

90 

d/B 

Input  Data  Rate 

0.602 

32 

M  Bits/s 

Size 

0.3 

0,  > 

19"  racV  *  6 1 

V 

■  .  .  .  ■■■  j 

1  Analyzer,  operating  withiut  external  »amy<e  .  no  held  system 

2  Maximum  continuous 


4.  4.  3  Fast  Fourier  Transform  Pi-rjys  :jn_n£ 

FFT  processing  to  compute  spectra)  mean  and  width  inay  be  performed 
either  by  computer  or  by  iipeciai  purpose  hardware.  Computer  processing  re¬ 
quires  that  the  weather  radar  signal  be  recorded  and  analyzed  in  non-real  time, 
but  subject  to  this  restriction  can  provide  essentially  any  desired  frequency  reso¬ 
lution.  The  PPP  is  compared  below  to  a  special  purpose  hardware  FFT  recently 
developed  by  Raytheon  for  pulse  compression. 


~>ft 

ppp 

Unit 

Range  Resoluth- . 

150 

i 

150 

Meter 6  (m) 

Rank  •>  elements 

2(1024} 

i 

1024 

- 

Velocity  Resolution 

0.5  >£  /204b 

i 

i 

0.5  >f  / 1024 
o 

m/  e 

Dynamic  Range 

i 

90 

dh 

Input  Data  Rate 

(321 

i 

i 

32 

M  Bit  s/ s 

Size 

2(4) 

i 

_ 1.  _ 

a.  3 

1 9"  rack  x  6  1  ] 
1—  . _J 

Numbers  in  (  )  indicate  anticipated  valuer  In r  a  1024  e,uti-  machine. 
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4. 5  Conclusions 


A  hardware  PPP  capable  of  computing  spectral  mean  and  width  for  each  of 
1024  1/2  -M-s  radar  range  gates  appears  completely  feasible  at  this  time.  The 
processor  can  be  fabricated  primarily  from  available,  designed  components  and 
is  expected  to  occupy  no  more  than  24"  of  19"  rack  space.  No  other  signal  pro¬ 
cessing  device  now  available  can  approach  the  dynamic  range  and  data  throughput 
capabilities  of  the  PPP,  except  possibly  the  FFT,  which  would  cost  approximately 
10  times  more  than  a  PPP  of  comparable  resolution. 


Table  7  gives  a  summary  of  the  design  goals  for  the  proposed  signal  pro¬ 


cessor. 


In  the  design  phase  of  the  PPP,  detailed  design  tradeoffs  of  the  following 
items  must  be  accomplished. 


1.  The  dynamic  AGC  loop 

2.  The  digital  function  generation 

3.  The  radar  calibration  system 

4.  Tnp  noi6e  monitoring  system 
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5.0  CONCLUSIONS  AND  RECOMMENDATIONS 


The  preceding  sections  of  this  report  show  analytically  and  verify  by  data 
from  controlled  simulation  experiments  that  the  pulse  pair  technique  produces 
estimates  of  spectral  mean  and  width  that  have  predictable  accuracy  over  a  wide 
range  of  non-ideal  operating  conditions,  Performance  for  unimodal  spectra  is 
insensitive  to  spectral  shape  and  the  degree  of  correlation  between  successive 
pulse  pairs;  the  theory  stands  up  under  conditions  of  poor  signal-to-noise  ratio. 

From  a  processor  standpoint,  it  is  shown  that  a  flexible  digital  pulse  pair 
processor  capable  of  handling  a  large  number  of  range  cells  simultaneously  is 
feasible, and  that  it  could  be  arranged  to  easily  match  a  variety  of  radars  and 
radar  operating  conditions,  including  PRFs,  pulse  widths  and  data  rates.  A 
theoretical  basis  for  optimising  the  waveform  parameters  of  a  radar  for  pulse 
pair  processing  is  also  presented. 

A  hardware  implementation  of  a  flexible  digital  pulse  pair  processor  suitable 
for  exploring  the  utility  of  the  technique  in  real  time  processing  of  weather  radar 
data  is  proposed.  The  design  goals  are  defined,  and  its  performance  is  compared 
with  other  contemporary  signal  processors. 

The  experimental  work  together  with  the  simulation  studies  demonstrate 
that  the  pulse  pair  technique  can  yield  good  quality  estimates  of  mean  and  variance 
usro?  rxdar  data  at  signal-to-noise  ratios  obtainable  with  existent  Doppler  weather 
radars,  f h-  n»;dware  tradeoff  studies  as  well  as  the  comparative  performance 
characteristics  demonstrate  that  the  pulse  pair  technique  is  not  only  feasible  but 
e! fie  t<*M. 

Bated  on  these  conclusions,  the  construction  of  a  flexible  experimental 
rnt>4fl  ot  such  a  real-time  pulse  pair  processor  is  recommended  for  the  purpose 
of  coliectii.g  a  large  data  base  on  observable  conditions.  Experience  gained  while 
conducting  research  with  the  experimental  model  will  provide  the  basis  for  defining 
tV,e  *prci£:caiions  for  simpler  models  intended  for  operational  use.  A  study  should 
also  be  made  to  determine  the  best  means  to  display,  plot,  and/cr  record  the  re¬ 
sulting  data  .or  further  analysis,  reduction  and  evaluation. 
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The  fabrication  phase  for  the  experimental  model  should  be  preceded  by  a 
design  evaluation  phase,  during  which  the  final  parameters  for  the  equipment  can 
be  established  from  the  results  of  tradeoff  studies. 

It  is  further  recommended  that  some  additional  analysis  be  performed  to 
provide  theoretical  guidelines  for  comparing  pulse  pair  analysis  to  spectra  data 
processors.  In  particular,  although  the  theory  of  pulBe  pair  estimation  has  now 
been  established,  no  theoretical  foundation  exists  at  preeent  for  predicting  the 
performance  of  width  estimates  derived  from  spectral  analyzer  data;  such  a 
theoretical  foundation  should  be  developed.  Analyses  are  also  required  to 
properly  relate  and  interface  the  pulse  pair  processor  with  particular  radars, 
ancillary  devices,  and  operational  usages  as  a  research  tool. 


5-2 


.ifMK^jrr-  ^  iiml|nT1|K»ui^r3l I f ir’C WJffi1, III1  ??'W!!3fflit 


REFERENCES 

1.  Hofstetter,  E.  M.  ,  "Simple  Estimates  of  Wake  Velocity  Parameters," 
Tech.  Note  1970-11,  23  April  1970,  M.  I.  T.,  Lincoln  Laboratory, 
Lexington,  Mass. 

2.  Rummler,  W.  D.  ,  "Introduction  of  a  New  Estimator  for  Velocity 
Spectral  Parameters,"  Tech.  Memo  MM-68-412  1-5,  3  April  1968, 

Bell  Telephone  Laboratories,  Whippany,  N.  J. 

3.  Rummler,  W.  D. ,  "Accuracy  of  Spectral  Width  Estimators  Using 
Pulse  Pair  Waveforms,"  Tech.  Memo  MM-68-412 1-14,  29  October 
1968,  Bell  Telephone  Laboratories,  Whippany,  N.  J. 

4.  Rummler,  W.  D.  ,  "Two-Pulse  Spectral  Measurements,  "  Tech.  Memo 
MM-68-4 12 1-15,  7  November  1968,  Bell  Telephone  Laboratories, 
Whippany,  N.  J. 

5.  "Doppler  Signal  Processing  and  Instrumentation  for  Modified  'Porcupine' 
C-Band  Pulse  Doppler  Radar,"  Final  Report,  May  1971,  on  Raytheon, 
Wayland,  Mass.  Contract  No.  F19628-68-0345  with  AFCRL  (pp  1-39 
through  1-49). 

6.  Berger,  T. ,  "Analysis  of  Data  from  Pulse  Pair  Feasibility  Experiment," 
Tech.  Memo  TB-108,  September  9,  1970,  Raytheon  Co.,  Wayland,  Mass. 

7.  Levin,  M.  J.,  "Generation  of  a  Sampled  Gaussian  Time  Series  Having 
a  Specified  Correlation  Function,"  IRE  Trans,  on  Information  Theory, 
Vol.  IT-6,  pp.  545-548,  December  I960. 


APPENDIX  A 


ANALYSIS  OF  ACCURACY  OF  SPECTRAL  PARAMETER  ESTIMATES 

BY  PULSE  PAIR  ESTIMATION 

H.  L.  Groginsky 


ABSTRACT 

This  appendix  gives  a  mathematical  analysis  of  the  performance  of  a 
pulse  pair  estimator  of  spectral  mean  and  variance.  It  assumes  that 
complex  data  (I  &  Q)  is  collected  at  a  uniform  rate  from  a  distributed 
target  whose  underlying  spectral  density  function  is  G(f)  and  that  G(f) 
is  observed  only  after  it  is  mixed  with  white  additive  noise. 

I.  Introduction 


If  the  signal  spectrum  is  written  as 


with  R(f)  real  and  positive  for  all  f. 


R(-f)  =  R(f), 

(2) 

f  R(f)  df  =  1. 

(3) 

/ f  R(f)  df  =  0, 

(4) 

y?  R(f)  df  =  i. 

(5) 

S  =  signal  power, 

and  if  the  data  is  given  by 

■k  =  *k  +  nk 

(6) 
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than  the  pulse  peir  estimators  of  f  and  W  are  given  by 


where  M-l 

X  '  M  ^  rk  rk+l 
n=0 

M-l 

Y  Irk!2 

n=0 

°a  *  <“k“k*>=  < !  ”kl  2  >  ' 


Z,  Spectral  Mean  Estimates 
Let 

oo 

g(r)  =  J  G(f)  ejZTrfT  df 
-00 

then 

<*k,k+i>= 

<’k,k+j>  ■  0 
<°k  nk+j  >  -  N  \s 
<"knktj>  ■  0 


where 

sk  =  s(kT). 

and  <},  .  is  the  Kronecker  delta  function, 
kj 


(7) 

(8) 

(9) 

(10) 

(ID 

(12) 

(13) 

(14) 

(15) 

(16) 


A-2 


Note  that 


J2nfT 


=  Se 


T  r  wT 


/.£*■) 


5  1 

1  ^ 


j  2wf  T 

Se  p  (wT) 


and  that  p  (wT)  is  real. 


M-l 

<x>  »  l  E 


rk  rk+l 


=  g(T) 


j  2iTf  T 

=  Se  °  p  (wT) 


Since  S  and  p  are  real. 


arg  <  X>  =  2  it  fQT 

which  shows  that  the  rrean  estimator  is  unbiased. 
Now  consider 


with  9  =  2tt  f  T 

o 

A  A 

If  6  is  to  be  a  useful  estimator,  8  *  8^,  where  8^,  is  the  true  mean 


phase  change  per  pair.  Writing 


£  =  eT  +  6§, 


Re  {X}  =  Re  {g(T)>  +  c  v 
Im  {X}  =  Im  {g(T)}  +  t2 


A-3 


1 


we  have 


lm{ g(T)}  +  e2 


cos  8t  Ro{g(T)}+ei 


Im  fa(T» 

-  ko  fiWj 

from  which  it  follows  that 


im7g(T)}j 


1  - 


Re  {g(T)}  e  -  Im  (g{T)}  e  , 

^  C  1 


since  cos  0T  =  . 

T  |  g(T)  | 


Now  by  defining  a  complex  error 


C1  +  jt2 


wo  may  write 


(T)  +  R*(T))<e-c*) 


4j  |  g(T)  j 


5*<T)  e  -  e*  g(T) 


2j  |  g(T)| 


* 

Im(g  e) 


=  Im 


ter) 


Now  define  a  new  complex  estimator  error 


£  1  \ 
Re  {g(T)y 


r*(T)l  U  -e  *) 
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Now  we  can  evaluate  the  requisite  expressions  using  the  assumption  of 
Gaussian  statistics  for  the  raw  input  data  as  follows: 


M-l  M-l 


«  ava  — *  in-*  ^  ^ 

1  X|  =  7?  k?0  So  <rkrk+1  ^+1> 


M-l  M-l 

r 


M  k=0 


£<  rk  rk+l  ><ri  ri+l  > 


+  <rkri  ><rk+l  ri+l> 


+  <  rk  Tm  ><'rk+l 


-:>] 


(39) 


From  (13)- (16)  it  follows  that 
term  1  *  |g(T)|2 

I  |  g  ( (i-k)T) 
term  2  =  j  , 

|  |  S  +  N  |  2 

term  3  *  0 . 

Note  s/N  is  the  per  sample  signal  to  noise  ratio.  For  not&tion&l  con¬ 
venience  we  will  write 

gn  =  g  (nT) .  (41) 


i  *  k 

i  =  k 


Then 


M-l  M-l 

<lxl*>a  I  I  2  +  "tE  £  I  gi.k|  2 

1  M  k=0  fco  1  K 


1  m-i  r  2  z] 

+  3T  E  l(S  +  N)  ' s  J 

M  k=0 

=  18^  2  ♦  TT  [<S  +  N)2-S2] 


M-l 


E 


where 

fM|n 

Similarly 

<x2> 


gJ  1  M 


(n) 


n=-(M-l) 


=  1  "M  “ 

IilL 

M 

,  M-l 

*  1J  E 

M  k=0 

M-l 

E 

£=0 

<rk 

M-l 

M-l 

s  -V  z 

E 

r<rk 

M  k=1 

£=l 

r 

*  * 
k  ‘ k+1  ri  T£+l 


(42) 


(43) 


+  <rk  ri>  <rk+l  r%+l> 


+  <rk  ri+l>  <r 


k+i  ri  >] 


(44) 
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Then  (13)  -  (16)  yields 
term  1  =  g2 


term  2=0 


term  3 


gi+l-kI*i*k“1  2»<  k-1.  k+1 


(S  +  N)  g_2 
g2  (S+N) 


l  =  k-1 


l  -  k+1 


and  substituting  these  in  (44)  yields 

M-l  M-l 

<X2>  =  g2  +  i-2  E  £  gi+1.k  g!k, 

1  M  k=0  i=0  X  k  1 


N(  g2  +  g*2]  (M-l) 


=  gl  +  M  £  8  n+l  gn-l  fM  (n) 

n=-(M-l) 


+  M  (  1  ’  M)  (g2  +  g.2) 


Now  substituting  (44)  and  (45)  in  (38  bis)  yields 


v«  {.«}  .A,  Rel^-s* 


N  (1-  -Jj) 


(8?  +  g  ,  ) 


♦  £ 

n=-(M-l) 


[  *n-l  \ 


•m"'1  j 
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For  large  M  and  small  wT  (46)  yields  approximately 


Var  {66} 


“  2M 


(1  4  N/S)2-! 
p(wT)|  2 


N_  2  p  (2tr) 

S  |  P  (wT) 


+  (2ttwT) 

|  P(wT)|  2 


1 

N 

S 

^•+2^1-  p  (2wT)  ^  +  (2tiwT) 

2M 

- 

|  P  (wT)  |  J 

(47) 


Equation  (47)  is  a  general  expression  for  predicting  the  variance  of  the 
pulse  pair  estimate  of  the  mean  Doppler  frequency  accounting  for  the  correla¬ 
tions  in  the  data  itself.  It  is  virtually  independent  of  actual  spectral  shape. 

3.  Width  Estimator 


We  have  shown  that  the  spectral  width  estimate  is  given  by 


A  2 
w 


[‘••V-Vl  • 


B 


(2  TT  T) 

If  the  estimate  is  useful,  it  must  have  small  errors  and  hence 
2 


w 


+  2  w  ( 6  w )  =  Bt  +  »?B 


(48) 


BT  + 


*lr  £1  +  gjlm  £ 

i  Sj  I  8 


gx  I  V 

~T~ 


(2itT  y 


where 

1*1 


A 


+  <i) 


+  <8lIm)  +  £2^ 


(50) 


and  Y  =  S  +  N  +  V 


(51) 
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Note  that 


«lr  h+  S1Irnt2  =  Re(g*fc) 


ec  that 


(2itT)2  W  (6w)  = 


It  follows  that 

\ 2 


[?  —  f <N]  ■ 


(52) 


(53) 


(2irT)4  ^  i^T")  Var  ^  =  Var  (6^)  +  Var(-~ -)-  2  cov  |(j^)  ( )|  (54) 


From  (32)  and  (33)  it  follows  that 
V.r<6*)=i  Re  |<|J-|>+ <-£_>! 


-  -I  (  jS+Nj 

"  2M  \  | 

i  I 


2-s2 


+  N 


1  v  («2  +  S.2) 


M-l 

+  £ 

n  =-(M-lp  |  g 


('in 
■lA  I  8, 


2 

T 


Following  the  methods  used  earlier  we  have 

M  M 


gn+l  gn- 


w  2  *  r  i  M  M 

<(S,>=  ?  [i?  £„  £,  <1  *1,1  *  I  *,l  *>-«♦««* 


(55) 


(56) 
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L[ 


n.-(M-l) 


r, „  i z  / 

f  1  .  2  \ 

,  8n+l  8n-l 

4  8n  Cl  1 

[UJ  ( 

M  \  \ 

s2"/ 

+  - 2 - 

81 

Sg!  J 

IM<n) 


This  is  the  requisite  expression  for  the  variance  of  the  width  parameter 
estimator. 

4.  Example 


To  reduce  the  equations  to  useful  form,  it  is  of  interest  to  evaluate 
them  for  a  Gauss  shaped  spectrum.  In  this  case 

1  -  *2 
1  e 


R(f)  = 


so  that 


p(t)  =  e 


M)Z 


S/w 


- 


g  =  S  e 
&n 


j2n  £  T 
o 


P  (nwT) 


j  2ir£  T 
o 


-  (2ir  nwT) 
2 


N  (g2  +  g  _2> 


,  N  !—ZJ 

S  |»|Z 


(63) 


®n±l  gn-l 


i  2 

PJ  * 


,  * 

gl  *  8.1 


=  2, 


gn  *„-l 


-  (2ttTw) 


.2  ♦  (n-l)2 


UttTw)'1'  n  (n-1) 


-  (2ttTw r  (  m+  1/2)  (m  -  1/2) 


1  (ZttTv,)2  -  (2TrmTw)2 
4 


which  leads  to 


2M  Var  {0}  = 


<  l  +  N/S)  -  1  -  2N./S  |  P  |  (1 


!"nl  2  <l-  I'll2)  'm1"'1 


n=-(M-l) 
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•  if  Hipty  ST1”-' 


2M  (2itT) 


'te) 


(1  +  N/S )2  -  1  1  +  2  IP 


M-l 

E  lPn|2[lt3lP. 


n=- ( M-l) 


+  2  T  lp2l  -  I  Pxl  2 


r  (2nTw)2  . 


These  expressions  may  be  simplified  as  follows  (for  large  M  and 
2xTw  <<1 


(  1  +  N/S)2  -  1  -  2  N/S  j  p2  |  (1  .  JL  )  -  n/s  (2  +  N  S)  -  2  N/S  |  p. 


=  N/S  [(N/S  +  2  (1  -  |  p-  |  )] 


2  ,Pn|2  (1-  lPl!2)fM(n) 

00 

«  I  Pni  2  < 1  -  I  pxI  2) 


w  (2irWT) 


»  2irwT 


=  2rrwx 


-(2TrnwT) 


•*>  2 

/•" 
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Similarly 


-][ 


+  n/sj  -  1  1  +  Z|p  I  *  +  |p2| 


8N,  ,2 

-T I  pil 


=  N/S  (N/S  +  2)  (  l  +  2|  p^|  2)  +  ]  p^|  -  ^gr-  |  p^|  2 


=  N/S  £  N/S  (l  +  2|Pl|2)  +  2  +  2 1  p2 1  -4|Pl|2j 

=  N/S  £  N/S  (l+2|Pl|2)  +  2(l-|Pl|2)2j  , 

„  r  „  I  (2^wT)Zl 

|pnl  [l  +  3lpl!  *4lpll  e  JfM(n) 

2  T  (2rrwT)2  -|(2™T)2 

=  V  Ip„I  1  +  3«  -4e  J 


2  T  (2rrwT)2  -4(2wT)2 

=  53  lf>nl  Ll  +3e  ~4e  -* 

«  |  (Z^twT)4  53  ipn'2 

«  |  (2ttwT)3  . 

As  a  re  suit  it  is  possible  to  reduce  the  variance  expressions  to 


2  M  Va 


r(S)  =  |  (  £♦*(».. -S 


+  2ttwT 


A- 15 


2M  (2nwT)2  Var  (2irwT) 


N 

r 


N  / 

*( 


1  +2  o 


■  (2irwT) 


^+2(l.e*(ZirwT)^ 


4-  ~  (2ttwT)3  ^tT 


(74) 


Figures  1  end  2  are  essentially  plots  of  (73)  and  (74)  and  enable  the 
prediction  of  system  performance  for  various  choices  of  parameters. 
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APPENDIX  B 

MATHEMATICAL  DERIVATIONS  RELATING  TO  THE 
GENERATION  OF  SIMULATED  DOPPLER  RETURN 


a.  Second  Moment: 


H(s )  H(.s)(-a*)  ds 


C  H 

2nj  *. 

-j* 


(a)  H(-a)  d 8 


j  j 

-l-  \  H(a)  H(-a)  da  -  — .  \  2  d: 

2ttj  J  .  2itj  J.  (s+a)  (■  s+a  f 


±  r _ i _ i 

da  (.a  +  a)  2  J 


a  =-  a 


Similarly, 


a  .  where  a  =  2ttw 


H(a )  H(-a)  (-a6)  ds  = 


i  a> 
4 


B-l 


i  .  -.ii'-.tt' .  . . . 


therefore  , 


=  w2 


b. 


Half  Power  Points  (w^): 


H(j„ 


a'' 


l 

TT  3  r 
J3  a  n/2 


w,=  a  -l  ] 


>/« 


Recursion  Relation: 


mi  -  H (v\ 
X!Z)  •  H(Z) 


7  *1  2  <-p 

Z  ale 


(l-2e  ^  z  1  *  e 


aT 


-2aTz‘l) 


Y(Z)  =  a*  Te  '  X(Z)  Z‘‘  +  2'  &T  Y(Z>  z’’  -  e'2aTY(Z)  Z-* 


which  yields 


a  _  -aT  -aT  -2a’T 

yn r  a  T  e  *n  i +  2  e  y„re  y 

n  n  - 1  n-l 


n-2 


d.  Signal  to  Noise: 

2  ^  yf 

-  Signal  Power  =  —  $  i  H(Z)  H(l/Z)  ~ 

2tt_j  J  Z 


The  factor  2  is  due  to  the  power  from  both  the  real  and  imaginary  channels. 
The  contour  of  integration  is  the  unit  circle.  The  Z-transiorm  is  used 
because  a  digital  signal  is  the  desired  output  from  the  simulator;  therefore, 
signal  power  should  be  found  in  terms  of  the  Z-transform.  To  continue, 

Z  dZ 


S  =  —  <§,  — 

2"j  J*  .at, 


Z -  e  *) 4  (Z-eaT)‘ 

Since  there  is  only  one  (second  order)  pole  inside  the  unit  circle  at  Z  =  e 
the  value  of  S  is  given  by: 


■aT 


f»  t  •  *44 ?sr3f  J  if/ »^^.T~-  ~.~: 


.  2nj  d  f  a4  T*  Z  "| 

2TTj  dZL(Z-eaT)2J 


Resulting  in 


=  Za  T‘ 


-aT  aT 
(3e  +  e  ) 

aT  - aT  , 
(e  -  e  )  J 


I  ' 


Noise  power.  N.  is  given  by: 


N  =  2  a 


The  factor  2  results  from  noise  power  in  the  real  and  imaginary  channels. 
Therefore  the  S/N  ratio  is  given  by: 


S  /  N  = 


A  '  -  aT  aT 

a4  T'  (  3^  ■-  e  ) 

2  -aT  , 

a N  <e  e  >5 


B-3/B-4 


'1*  Mflttl'N1  ■"  *  '">?!■  ■■■»  *-  •  nr 


APPENDIX  C 


FORTRAN  MAIN  PROGRAMS 


C-l 

DOPGEN 

C-2 

PPSTAT 

C-3 

SPEC 

COMMENT — SYSTEM  SUPPLIED  SUBROUTINES  RELUIREG —  DQPOLT  >  ENDPLT 
COMMENT — USER  SUPPLIED  SUBROUTINES  REQUIRED — G'.  RN»  FORT 
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PROGRAM 


TRACE 


COC  6600  ETN  V  3. 0-P27  0 ♦  0PT*0  03/2 


PROGRAM  SPEC (INPUT, OUTPUT, TAPE5=INPUT,TAPE6=OUTPUT,TAPE8,TAPE9> 


COMM; NT —  SYSTEM  SUP»LIE0  SUBROUTINES  REQUIRED —  OQPLOT,ENOPLT 


FTN(A) 


SETCORE. 


REQUEST, TAPE33. 


REQUEST »TAPE3,HI.(RPK08/NORING> 
i  r.n. 


UNLOAO,TAPE3}. 


UNLOAD, TAPES. 


REAL  MEAN 


DIMENSION  Z (10241# A (1024)  , X  ( 10 24 > , Y ( 1S3&) 


N=1024 


I  4JBtl*?il 


I °LOT«0 


I K  '•■■Mil***# 


10  JJUUU I  ia-LXMlIi,  HOC,  I  PL  OX,  MI  AS  . 

to,  1000)  IDlX,  NUNJT  ,1.  IM1T1  ,L  IMIT2,IRN  ,  IPlOT,  I  BIAS 

-  .  ..  j:r.U8uc^ev(UGa  ..ta  .jlu  _ _ 

I'  » LIMIT  \.Vi.  3)GQ  Tv->  \i 

- .QjQ-J.U.^JL.. - - - 

i?  uIMITI-*: 

.  ..  ..•.IMlT.aAl6.2V-  - . —  _ _ _  _ _ 

13  CONTINUE 


Irf  IRM.ET.1MEK* 

15  PT*lMHUN.r.  —  «  .Ft  3  ,r  1 1,12,  13,  (Z(U  ,1*1,1024) 

. -  -  -  .  — . - 

:■-<  Tc-T.N',  To 

O'  ?0  I  ■*  1  ,  t  C  L‘  4 

■  Z-t-ilLAZl.Ll.YAU-l _ _ _ _  _  .  _ 


PROGRAM 


SPEC 


TRACE 


roc  6600  FT N  V3.0-P270*  OPT»0  03/2 


CALL  F0RT(Z,1, S, IFS ,IFERR> 


00  30  1*1,1024 

JO.  Y(H=REAL  <Z(I)  )**2»AIMAG(Z(I))**2 

_ YMAX=Y(1) _ 

00  40  1=1,1024 

40  IF<THAX.LT.Ym>  YHAX»Y<I) _ 


_ IF{IPt.OT,£Q,a)  60  TO  75 

00  50  1=1, 102* 

50  rcn=Y(i>  /ymax _ 


_ fff*-i,q _ 

00  60  1=1,1024 

_ FFF  =  FFF»i.q  _ ■  ■■  ■  _ _ _ 

X (I)=FFr 

60  Y  (I>=lQ,»ALQSlQ(Y(in _ 

- CALL.  0Q?LQT..(LA.a.« X«  Y  ,N ,  1.H  .1 , 12M2 0 . LQ G  (HAS)  .12,0,  m 

K PLOT  =1 


00  70  J= 1 , 10  24 

7-0.  .-LU).=  10« * T.1Y..U.1./10, J».YMAX _ I 

75  CONTINUE 


I F  ( IA  VG,  EQ.  01,0  TO  20  0 

_  _ 

JO  76  1=741,940 

76-AVC;:AVQ*YlIl _ 

AYG=AVG/100. 


DO  77  1=1,1024 

--gz_Y.(i.>=Ym’4ys _ 

200  I F ( IRIAS*  EC, J > GO  TO  230 

_ LSrt£I-512-IBI4S _ 

00  210  1  =  1,1024 

_ J1=1924»ISHFT»1-I _ 

Kl=1024*l-I 

218  YIJ1)=Y(<1) _ 

00  220  I  =  l,I3-tFT 

_ J2=ISHFT  tl*I _ 

K 2=1024*1 SHFT  *1-1 

-22Q  Y  1 J 2 ?  =Y(K2) _ 

230  CONTINUE 

_ WRITE  16. 20001 _ 


00  90  JsLIHITl,LlMIT2 
FH=FM*1. _ 


IF(r<JI.'.E.YHIMV»30  TO  90 


Y  SQg.Y.t  J.1 _ 

S9=S3*YSQMF1-UKP>  *•? 
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